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COMPLEX  ENVELOPE  PROPERTIES,  INTERPRETATION, 

FILTERING,  AND  EVALUATION 

INTRODUCTION 

When  a  narrowband  input  excites  a  passband  filter,  the  output 
time  waveform  y(t)  is  a  narrowband  process  with  low-frequency 
amplitude-  and/or  phase-modulations.  The  evaluation  of  this 
output  process  y{t)  can  entail  an  extreme  amount  of  calculations 
if  the  detailed  behavior  of  the  higher-frequency  carrier  is 
tracked.  A  much  better  procedure  in  this  case  is  to  concentrate 
instead  on  determination  of  the  low-frequency  complex  envelope  of 
the  narrowband  output  process  y(t)  and  to  state  the  carrier 
frequency  associated  with  it.  Then,  the  detailed  nature  of  the 
output  can  be  found  at  any  time  points  of  interest  if  desired, 
although,  often,  the  complex  envelope  itself  is  the  quantity  of 
interest . 

The  complex  envelope  of  output  y(t)  is  determined  from  its 
spectrum  (Fourier  transform)  Y(f)  by  suppressing  the  negative 
frequencies,  down-shifting  by  the  carrier  frequency,  and  Fourier 
transforming  back  into  the  time  domain.  For  a  complicated  input 
spectrum  and/or  filter  transfer  function  with  slowly  decaying 
spectral  skirts,  these  calculations  can  encounter  a  large  number 
of  data  points  and  require  large-size  fast  Fourier  transforms 
( FFTs )  for  their  direct  realization.  In  this  case,  the  use  of 
collapsing  or  pre-aliasing  [1;  pages  4-5]  can  be  fruitfully 
employed,  thereby  keeping  storage  and  FFT  sizes  small,  without 
any  loss  of  accuracy.  This  procedure  will  be  employed  here. 
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As  will  be  seen,  when  the  complex  envelope  is  re-applied  to 
the  one-sided  carrier  term  and  the  real  part  taken,  the  exact 
narrowband  waveform  y(t)  is  recovered.  However,  if  the  complex 
envelope  itself  is  the  quantity  of  interest,  it  has  some 
undesirable  features.  The  first  problem  is  related  to  the  fact 
that  if  waveform  y(t)  has  any  discontinuities  it  time,  its 
Hilbert  transform  contains  logarithmic  infinities,  which  show  up 
in  the  complex  envelope.  The  second  problem  is  generated  by  the 
operation  of  truncating  the  negative  frequencies  in  spectrum 
Y ( f ) ;  this  creates  a  discontinuous  spectrum  which  leads  to  a  very 
slow  decay  in  time  of  the  magnitude  of  the  complex  envelope. 

Since  numerical  calculation  of  the  complex  envelope  is 
necessarily  accomplished  by  sampling  spectrum  Y(f)  in  frequency  f 
and  performing  FFTs ,  this  slow  time  decay  leads  to  significant 
aliasing  and  distortion  in  the  time  domain  of  the  computed 
quantities . 

Because  these  features  in  the  mathematically  defined  complex 
envelope  are  very  undesirable,  there  is  a  need  to  define  and 
investigate  a  modified  complex  envelope  which  more  nearly 
corresponds  to  physical  interpretation  and  utility.  The  time 
discontinuities  in  y(t)  show  up  in  Y(f)  as  a  1/f  decay  for  large 
frequencies,  whereas  the  truncation  of  the  negative  frequencies 
of  Y(f)  shows  up  as  a  discontinuity  directly  in  f.  Both  of  these 
spectral  properties  can  be  controlled  by  filtering  the  truncated 
spectral  quantity,  prior  to  transforming  back  to  the  time  domain, 
'.'e  will  address  this  filtered  complex  envelope  and  its  efficient 
evaluation  in  this  report. 
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When  the  waveform  y(t)  is  real  and/or  causal,  its  spectrum 
Y ( f )  possesses  special  properties  which  enable  alternative 
methods  of  calculation.  Thus,  it  sometimes  suffices  to  have  only 
the  real  (or  imaginary)  part  of  Y(f)  and  to  employ  a  cosine  (or 
sine)  transform,  rather  than  a  complex  exponential  transform. 

The  aliasing  properties  of  these  special  transforms,  when 
implemented  by  means  of  FFTs,  will  also  be  addressed  here. 
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ANALYTIC  WAVEFORM  AND  COMPLEX  ENVELOPE 

Waveform  y(t)  is  real  with  amplitude  modulation  a(t)  and 
phase  modulation  p(t)  applied  on  given  carrier  frequency  f  ; 
however,  y(t)  need  not  be  narrowband.  That  is, 

y(t)  =  a(t)  cos[2nfQt  +  p( t ) ]  =  Re{z(t)  exp(i2nfQt)l  ,  (1) 

where  complex  lowpass  waveform 

z(t)  =  a(t)  exp[ip(t) ]  (2) 

will  be  called  the  imposed  modulation.  The  corresponding 
spectrum  of  imposed  modulation  z(t)  is 

Z(f)  =  |  dt  exp(-i2nft)  z(t)  .  (3) 

(Integrals  without  limits  are  from  -®°  to  +».  )  The  magnitude 
of  spectrum  Z(f)  is  depicted  in  figure  1;  it  is  generally  concen¬ 
trated  near  frequency  f  =  0.  The  graininess  of  the  curves  here 
is  due  to  plotter  quantization,  not  function  discontinuities. 

From  (1),  since  waveform 

y(t)  =  \  z ( t )  exp(i2nfQt)  +  ^  z*(t)  exp(-i2nfQt)  ,  (4) 

its  spectrum  can  be  expressed  as  (see  figure  1) 

Y(f)  =  \  Z ( f  —  f o )  +  \  Z*(-f-fQ)  ;  Y ( - f )  =  Y*(f)  .  (5) 

It  will  be  assumed  here  that  y(t)  has  no  dc  component;  that  is, 
Yff)  contains  no  impulse  at  f  =  0. 
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Figure  1.  Spectral  Quantities 
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ANALYTIC  WAVEFORM 

The  single-sided  (positive)  frequency  spectrum  is  defined  as 

Y+(f)  =  2  U(f)  Y(f)  =  U(f)  Z( f — f 0 )  +  U ( f )  Z* ( -f-f q )  =  (6) 

=  Z ( f-f o )  -  U( -f )  Z ( f  —  f o )  +  U ( f )  Z* ( —  f  —  f 0 )  for  all  f  .  (7) 

Here,  U(x)  is  the  unit  step  function;  that  is,  U(x)  is  1  for 
x  >  0  and  U(x)  is  0  for  x  <  0.  The  analytic  waveform 
corresponding  to  y ( t )  is  then  given  by  Fourier  transform 

Y+(t)  =  J  df  exp( i2nft )  Y+(f)  .  (8) 

In  order  to  further  develop  (8),  we  define  a  single-sided 
(negative)  frequency  function 

0  for  f  >  0' 

N(f)  =  U(-f)  Z ( f-f  )  =  ,  (9) 

Z(f-f  )  for  f  <  0 
l  o' 

wRich  can  be  determined  directly  from  the  spectrum  Z(f)  of  the 
imposed  modulation  z(t)  in  (2)  if  f  is  known.  The  magnitude  of 
N ( f ) ,  scaled  to  peak  value  1,  is  sketched  in  figure  1;  it  is 
small  if  f  is  large,  and  is  peaked  near  f  =  0.  The  complex  time 
function  corresponding  to  (negative  frequency)  function  N(f)  is 

0 

n(t)  =  |  df  exp(i2nft)  N(f)  =  J  df  exp(i2nft)  Z(f-fQ)  .  (10) 

—  00 

With  the  help  of  (9)  and  (10),  the  single-sided  spectrum 
Y+(f)  in  (7)  can  now  be  expressed  as 


7 


TR  8827 


Y+(f)  =  Z(f-fQ)  "  N<f>  +  N*(-f)  ,  (11) 

with  corresponding  analytic  waveform  (8) 

y+(t)  =  exp(  i2rtfQt )  z(t)  -  n(t)  +  n*(t)  »  (12) 

=  exp(i2nfQt)  z(t)  -  i  2  Im(n(t)j  .  (13) 

That  is,  the  analytic  waveform  is  composed  of  two  parts,  the 
first  of  which  is  what  we  would  typically  desire,  namely  the 
imposed  modulation  (2)  shifted  up  in  frequency  by  f  .  The  second 
term  in  (13),  which  is  totally  imaginary,  is  usually  undesired; 
it  can  be  seen  from  (10)  and  |N(f)|  in  figure  1  to  be  generally 
rather  small.  There  also  follows  immediately,  from  (13)  and  (2), 
the  expected  result 

Re I y+ ( t ) }  =  a ( t )  cos [ 2nfQt  +  p(t)]  =  y(t)  .  (14) 

Since  analytic  waveform  y+(t)  can  also  be  expressed  as 
y+(t)  =  y(t)  +  i  yH(t)  =  y ( t)  +  i  y(t)  ©  yy  =  y(t)  +  ijdu 

(15) 

where  yH(t)  is  the  Hilbert  transform  of  y(t)  and  ©  denotes 
convolution,  (13)  and  (2)  yield 

yH(t)  =  a ( t )  sin[2nfQt  +  p ( t ) ]  -  2  Im(n(t)|  .  (16) 

I f  we  define  (real)  error  waveform  e(t)  as  the  difference  between 
the  Hilbert  transform  of  (1)  and  the  quadrature  version  of 
original  waveform  (1),  we  have 
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e(t)  s  y H ( t )  -  a(t)  sin[2nfQt  +  p(t)]  = 

=  -  2  Im{n{ t ) }  =  i  [n( t)  -  n*(t) ]  = 

0 

=  -  2  Im  J  df  exp(i2nft)  Z(f-fQ)  = 

—  00 

=  -  2  Im  j  exp(  i2rtf  Qt )  J  df  exp(i2nft)  Z  (  f )  I  , 
'  —00  ' 


(17) 

(18) 

(19) 


(20) 


where  we  used  (16)  and  (10).  The  error  spectrum  is,  from  (18) 
and  ( 9 ) , 

E ( f )  =  i  [ N( f )  -  N* ( -f ) ]  =  (21) 


'-i 

1  2(£-fo» 


for  f  >  0 
for  f  <  0 


(22) 


Then,  E ( -f )  -  E*(f).  The  magnitude  of  E(f)  is  displayed  in 
figure  2;  it  is  generally  small  and  centered  about  f  =  0. 
The  total  energy  in  real  error  waveform  e(t)  is 


J  dt  [e( t) ]2  =  |  df  |E(f) |2  = 


=  J  df  [ N( f )  -  N* ( -f ) ]  [n* ( f )  -  N ( —  f ) ]  = 

=  J  df  [ | N ( f ) | 2  +  |N( -f ) | 2 ]  =  2  J  df  | N ( f ) | 2  = 

0  -fQ 

=  2  |  df  | Z( f-fQ)  | 2  =  2  |  df  |  Z ( f  )  |  2  ,  (23) 

—  00  —00 
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Figure  2.  Error  and  Complex  Envelope  Spectra 
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where  we  used  (21),  the  single-sided  behavior  of  N(f),  and  (9). 
This  is  just  twice  the  energy  in  the  spectrum  Z(f)  of  imposed 
modulation  z(t)  below  frequency  -f  ;  inspection  of  figure  1 
reveals  that  this  quantity  will  usually  be  small. 


COMPLEX  ENVELOPE 

The  complex  envelope  ^(t)  waveform  y(t)  is  the  frequency 
down-shifted  version  of  analytic  waveform  y+(t): 

Y(t)  ■  y+(t)  exp( -i2nfQt )  =  (24) 

=  z(t)  +  i  e(t)  exp(-i2nfQt)  ,  (25) 

where  we  used  (13)  and  (18)  and  chose  to  downshift  by  fQ  Hertz, 
the  known  carrier  frequency  in  (1).  Waveforms  z(t)  and  e(t)  are 
lowpass,  as  may  be  verified  from  their  spectra  in  figures  1  and 
2.  The  spectrum  of  the  complex  envelope  is,  from  (25), 

Y(f)  =  Z ( f )  +  i  E(f+f0)  .  (26) 

Equations  (25)  and  (26)  show  that  the  complex  envelope  and  its 
spectrum  are  each  composed  of  a  desired  component  and  an  error 
term. 

The  magnitudes  of  the  complex  envelope  spectrum  Y ( f )  and  its 
error  component  are  displayed  in  figure  2;  Y(f)  is  discontinuous 
at  f  =  -f  but  has  zero  slope  as  f  ■>  -fQ,  w^et^er  fr°m  above  or 
below.  The  left  tail  of  Z(f)  and  shifted  error  spectrum, 
i  E(f+fQ),  interact  so  as  to  yield  Y(f)  =  0  for  f  <  -fQ;  this  is 
most  easily  seen  from  a  combination  of  (24)  and  (6),  namely 
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Y(f)  =  Y+(f+f0)  =  2  U(f+fQ)  Y(f+fQ)  =  (27) 

=  U(f+fQ)  [2(f)  +  Z* ( -f-2f Q ) ]  .  (28) 

The  results  for  the  error  spectrum  and  energy  in  (22)  and 
(23),  respectively,  were  originally  derived  by  Nuttall  [2]; 
however,  we  have  augmented  those  results  here,  to  give  detailed 
expressions  for  the  error  and  complex  envelope  waveforms  and 
spectra.  There  are  no  approximations  in  any  of  the  above 
relations;  they  apply  to  waveforms  with  arbitrary  spectra, 
whether  carrier  frequency  f  is  large  or  small. 

EXTRACTED  AMPLITUDE  AND  PHASE  MODULATIONS 

It  is  important  and  useful  to  also  make  the  following 
observations  relative  to  the  amplitude  and  phase  modulations  that 
can  be  extracted  from  the  complex  envelope  x(t).  Define 

A ( t )  =  | y( t ) |  ,  P(t)  =  arg{y(t) }  .  (29) 

Then,  from  (14)  and  (24),  the  original  waveform  can  be  expressed 
in  terms  of  these  extracted  amplitude  and  phase  modulations  as 

y(t)  =  Re[y(t)  exp(i2nfQt)|  =  A(t)  cos[2nfQt  +  P(t)]  .  (30) 

However,  complex-envelope  modulations  A(t)  and  P(t)  in  (29)  and 
(30)  are  not  generally  equal  to  imposed  modulations  a(t)  and  p ( t ) 
in  (1),  as  may  be  seen  by  reference  to  (25).  Namely,  complex 
envelope  y(t)  is  equal  to  complex  lowpass  waveform  z(t)  in  (2) 
only  if  error  e(t)  is  zero.  But  the  energy  in  waveform  e(t),  as 
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given  by  (23),  is  zero  only  if  imposed  spectrum  Z ( f )  in  (3)  is 

zero  for  f  <  -f  Q .  When  Z(f)  is  not  zero  for  f  <  -f  ,  complex- 

envelope  modulations  A(t)  and  P(t)  do  not  agree  with  imposed 

modulations  a(t)  and  p(t),  despite  the  ability  to  write  y(t)  in 

the  two  similar  real  forms  (1)  and  (30)  involving  an  amplitude- 

and  phase-modulated  cosine  with  the  same  f  . 

o 

Another  interesting  property  of  form  (30)  is  that  its 
quadrature  version  is  identically  the  Hilbert  transform  of  y(t). 
This  is  in  contrast  with  the  quadrature  version  of  (1)  involving 
imposed  modulations  a(t)  and  p(t);  see  (17)  -  (20).  To  prove 
this  claim,  observe  that  the  quadrature  version  of  the  last  term 
of  (30)  is,  using  (29), 

q(t)  —  A(t)  sin[ 2nf Qt  +  P(t)]  =  (31) 

=  [ A( t )  exp(iP(t)  +  i2nf  t)  -  A(t)  exp(-iP(t)  -  i2nf  t)]  = 

1"  o  o 

=  jjtY(t)  exp( i2nfQt )  -  y*(t)  exp( -i2nfQt ) ]  .  (32) 

The  spectrum  of  this  waveform  is 

Q(f)  =  3^[v(f-f0)  -  Y*(-f-f0)J  =  -i[U(f)  Y(f)  -  U(-f)  Y*( -f ) ]  = 
(-i  Y ( f )  for  f  >  0] 

=  -i  sgn(f)  Y ( f )  =  Y  (f)  ,  (33) 

{  i  Y ( f )  for  f  <  Oj  H 

where  we  used  (27),  the  conjugate  symmetry  of  Y(f),  sgn(x)  =  +1 
for  x  >  0  and  -1  for  x  <  0,  and  (6)  in  the  form 

Y+(f)  =  2  U(f)  Y ( f )  =  [1  +  sgn(f)]  Y ( f )  =  Y ( f )  +  i  YH(f)  ,  (34) 
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the  latter  result  following  from  (15).  Thus,  (33)  and  (31)  yield 
the  desired  result 

YH(t)  =  ttE  ®  y(t)  =  q(t)  =  A(t)  sin[2nf0t  +  P(t)]  •  (35) 

This  simple  connection  between  (35)  and  (30)  holds  in  general 

when  modulations  A(t)  and  P(t)  are  extracted  from  the  complex 

envelope  according  to  (29);  there  are  no  narrowband  assumptions 

required.  The  more  complicated  connection  between  (16)  and  (1), 

which  is  applicable  for  the  imposed  modulations,  involves  an 

error  term;  this  error  is  zero  if  and  only  if  spectrum  Z(f)  in 

(3)  is  zero  for  f  <  -f  . 

o 

SPECTRUM  Y ( f )  GIVEN 

All  of  the  above  results  have  presumed  that  waveform  y(t)  in 
the  form  (1)  was  available  as  the  starting  point.  But  there  are 
many  problems  of  interest  where  spectrum  Y(f)  is  the  initial 
available  quantity,  rather  than  y(t).  For  example,  the  output 
spectrum  Y(f)  of  a  linear  filter  L(f)  subject  to  input  spectrum 
X ( f )  is  given  by  Y(f)  =  L(f)  X(f)  and  can  often  be  easily  and 
directly  computed.  In  this  case,  there  are  no  given  amplitude 
and  phase  modulations  a(t)  and  p(t)  as  in  (1);  in  fact,  there  is 
not  even  an  obvious  or  unique  center  frequency  for  a  given 
spectrum  Y(f).  Nevertheless,  many,  but  not  all,  of  the  relations 
above  hold  true  under  appropriate  definitions  of  the  various 
terms . 
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Given  spectrum  Y(f)  with  conjugate  symmetry,  Y ( —  f )  =  Y*(f), 
we  begin  with  its  corresponding  real  waveform 

y(t)  =  J  df  exp(i2rcft)  Y(f)  .  (36) 

The  Hilbert  transform  of  y(t)  and  its  spectrum  are  given  by 

yH{t)  =  FE  ®  y(t)  '  YH(f)  =  _i  s9n( f )  y(f)  •  (37) 

The  single-sided  spectrum  and  analytic  waveform  are,  respectively 

Y+(f)  =  2  U ( f )  Y ( f )  =  [1  +  sgn(f)]  Y ( f )  =  Y(f)  +  i  YH(f)  ,  (38) 

CO 

y+(t )  =  2  J  df  exp{ i2nft )  Y(f)  =  y(t)  +  i  yH ( t )  .  (39) 

0 

Up  to  this  point,  all  the  functions  are  unique  and  nothing 
has  changed.  However,  we  now  have  to  choose  a  "center  frequency" 
fc  of  Y+(f),  since  none  has  been  specified;  this  (somewhat 
arbitrary)  selection  process  of  fc  is  addressed  in  appendix  A,  to 
which  the  reader  is  referred  at  this  point.  Hence,  we  take  fc  as 
given  and  define  lowpass  spectrum 

Y(f)  =  y+(f+fc)  =  2  U(f+fc)  Y(f+fc)  .  (40) 

The  corresponding  complex  envelope  is 

¥(t)  =  y+(t)  exp(-i2nfct)  .  (41) 

We  define  the  complex-envelope  amplitude  and  phase 
modulations  as  in  (29): 

A(t)  =  | y( t) |  ,  P(t)  =  argly(t))  =  arg{y+(t)j  -  2nfct  .  (42) 


15 


TR  8827 


Then,  from  (39),  (41),  and  (42),  we  have 

y(t)  =  Re{y  ( t ) }  =  Re{y(t)  exp(i2nf  t) )  =  A(t)  cos[2nf  t  +  P(t)]. 

(43) 

Now  when  we  define  the  quadrature  version  of  the  right-hand  side 
of  (43)  in  a  manner  similar  to  (31),  but  now  employing  f  instead 
of  (unspecified)  f  ,  the  same  type  of  manipulations  as  in 
(31)  -  (35)  yield  relations  identical  to  those  given  above: 

Q ( f )  =  YH( f )  ,  yH(t)  =  q ( t )  =  A(t)  sin[2nfct  +  P(t)J  .  (44) 

Because  the  choice  of  center  frequency  f  of  single-sided 
spectrum  Y+(f)  is  somewhat  arbitrary  (see  appendix  A),  this  makes 
complex  envelope  y(t)  and  its  extracted  phase  P(t)  somewhat 
arbitrary.  However,  the  argument,  2itfct  +  P(t)  =  arg{y+(t)l,  of 
(43)  and  (44)  is  not  arbitrary,  as  seen  directly  from  (41)  and 
the  uniqueness  of  y+(t)  in  (39).  Furthermore,  extracted 
amplitude  modulation  A(t)  in  (42)  has  no  arbitrariness  since  it 
is  given  alternatively  by  |y+(t)|,  according  to  (41). 

Since  A(t)  and  P(t)  are  lowpass  functions,  we  can  compute 
them  at  relatively  coarse  increments  in  time  t.  Then,  if  we  want 
to  observe  the  fine  detail  of  y(t),  as  given  by  (43),  we  can 
interpolate  between  these  values  of  A(t)  and  P(t)  and  then 
compute  the  cosine  in  (43)  at  whatever  t  values  are  of  interest. 
This  practical  numerical  approach  will  reduce  the  number  of 
computations  of  A(t)  and  P(t)  required;  in  fact,  in  many 
applications,  A(t)  and  P(t)  will  themselves  be  the  desired  output 
quantities  of  interest,  rather  than  narrowband  waveform  y(t)  with 
all  its  unimportant  high-frequency  detail. 
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EXAMPLE 

Consider  the  fundamental  building  block  of  systems  with 
rational  transfer  functions,  namely 


y(t)  =  U ( t )  exp(-at)  cos(2nfQt  +  $ )  ,  a  >  0  ,  fQ  >  0  ,  (45 

where  U(t)  is  the  unit  step  in  time  t.  Let 


0)  =  2nf  ,  w  =  2nf  ,  c  =  a  -  io) 

'  o  o  '  o 


(46 


Then,  from  (1)  and  (2),  the  imposed  modulations  are 


a  ( t )  =  U(t)  exp(-at),  p(t)  =  4.,  z(t)  =  U(t)  exp(i^-at)  ,  (47 


yielding,  upon  use  of  (3)  and  (46),  spectrum 


Z(f)  =  =  gxP(^) 

a  +  i2nf  a  +  io) 


(48 


From  (5)  and  (48),  the  spectrum  of  y(t)  is 


v(f)  -  f 


exp(  i<)> ) 


exp(-i4>) 


a  +  i(w  -  wQ)  a  +  i(o)  +  o)Q)J  ' 


49) 


and  (6)  yields  single-sided  spectrum 


Y  (f)  =  u  (  f )  f - -  +  - exP(-^) 

+  l  '  1  ;(a  +  i(w  -  w  )  a  +  i(u)  +  o)  ) 


50 


Now  we  use  (9),  (48),  and  (4o)  to  obtain  (negative)  spectrum 


N  (  f )  =  U(-f)  = 

'  '  v  '  C  +  10) 


0  for  f  >  0 

SffiUil  for  £  <  0 
C  +  10) 


(51 
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Then,  from  (10),  the  corresponding  complex  time  waveform  is 
0  c 

n(t)  -  J  df  exp(i„t)  J  ^  exp(tu)  .  (52) 

-00  C-i00 

For  t  <  0,  let  x  =  |t|u  =  -tu,  to  get 

c  t| 

n(t)  =  S3£BU|zetl  ^  exp(-x)  =  ^  exp(i^-ct)  E^cltl)  ,  (53) 

c  1 1  -i°° 

where  E^z)  is  the  exponential  integral  [3;  5.1.1].  It  is 
important  to  observe  and  use  the  fact  that  the  path  of 
integration  in  the  complex  x-plane  in  (53)  remains  in  the  fourth 
quadrant  and  never  crosses  the  negative  real  x-axis 
[ 3 ;  under  5.1.6]. 

Also,  for  t  >  0,  let  x  =  -tu  in  (52),  to  get 

-ct 

n(t)  =  eXPti2nCt>  exP('x)  =  2^  exp(i*-ct)  Ej(-ct)  .  (54) 

-ct+i«> 

Here,  the  contour  of  integration  remains  in  the  second  quadrant 
of  the  complex  x-plane  and  again  does  not  cross  the  negative  real 
x-axis  f3;  under  5.1.6].  The  combination  of  (53)  and  (54)  now 
yields  complex  time  waveform 

n(t)  =  — ~  exp(i^-ct)  E^(-ct)  for  all  t  ^  0  .  (55) 
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Now,  we  use  (18)  to  obtain  real  error  waveform 

e(t)  =  -  ^  Re{exp(i<|>-ct)  E1(-ct)J  for  all  t  *  0  .  (56) 

(Or  we  could  directly  use  (20)  with  (48).)  The  corresponding 
error  spectrum  follows  from  (22),  (48),  and  (46)  as 

for  f  >  0 

for  f  <  0 

From  (16),  (17),  (47),  and  (56),  the  Hilbert  transform  of 
y ( t )  is 

yH(t)  =  U(t)  exp  (-at)  sin(2nfot+4>)  - 
-  ^  Re { exp( i$-ct )  E1(-ct)}  for  all  t  *  0.  (58) 

In  addition,  using  (15)  and  (45),  the  analytic  waveform  is 

y+(t)  =  U(t)  exp(i<fr-ct)  -  i  -  Re{exp( i$-ct)  E1(-ct)} 

for  all  t  *  0.  (59) 

The  complex  envelope  follows  from  (25),  (47),  and  (56)  as 

y(t)  =  U(t)  exp(  i<J>-at )  -  ~  exp(-iwQt)  Re  { exp(  i^-ct )  E1(-ct)} 

for  all  t  *  0  .  (60) 

The  corresponding  spectrum  is,  from  (27)  and  (50), 
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The  extracted  amplitude  and  phase  modulations  A(t)  and  P(t)  of 
complex  envelope  ^(t)  are  now  available  by  applying  (42)  to 
(60).  Since  the  first  term,  by  itself,  in  (60)  has  the  imposed 
amplitude  and  phase  modulations  a(t)  and  p(t)  as  specified  in 
(47),  A(t)  cannot  possibly  equal  a(t),  nor  can  P(t)  equal  p(t). 
This  example  is  an  illustration  of  the  general  property  stated  in 
the  sequel  to  (30).  The  reason  is  that  spectrum  Z(f)  in  (48)  is 
obviously  nonzero  for  f  <  -fQ. 

From  (23)  and  (48),  the  energy  in  error  waveform  e(t)  is 


—  00 


—  arctan 

n 

a 

(62) 


For  comparison,  the  energy  in  desired  component  z(t)  in  complex 
envelope  ^(t)  of  (25)  ^s'  from  (47), 


|  dt  |z(t) I2  =  2H  •  (63) 

—  00 
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SINGULAR  BEHAVIOR 

Since  [3;  5.1.11  and  footnote  on  page  228] 

E^(z)  =  -  ln(z)  -  y  +  Ein(z)  ,  (64) 

where  Ein(z)  is  entire,  the  error  waveform  in  (56)  has  a 
component 

-  i  Re { -exp( i$-ct)  ln(-ct)}  = 

=  ^  Re[exp(i<)>-ct)  [ln(-c  sgn(t))  +  ln|t|]  for  t  *  0  ,  (65) 

of  which  the  most  singular  component  is 

^  In  1 1 1  exp(-at)  cos(ooot  +  <t>)  ~  ^  cos(^)  In  1 1 1  as  t  ->  0  .  (66) 

The  only  situation  for  which  this  logarithmic  singularity  does 
not  contribute  an  infinity  as  t  ^  0  is  when  <}>  =  -  n/2  (or  n/2). 
That  corresponds  to  the  special  case  in  (45)  of 

y(t)  =  U ( t )  exp(-at)  sin(u)Qt)  for  <}>  =  -  n/2  ,  (67) 

which  is  zero  at  t  =  0;  that  is,  y(t)  is  continuous  for  all  t. 
However,  even  for  <}»  *  -  n/2  in  the  first  term  of  (66),  the 
product  ln|t|  sin(ooot)  has  an  infinite  slope  at  its  zero  at 
t  =  0,  leading  possibly  to  numerical  difficulties. 

The  spectrum  Y(f)  follows  from  (49)  as 

w 

Y  (  f  )  =  ~ 2 - 2  - 2  for  ♦  =  ~  §  ,  (68) 

a  +  w  +  i2aw  -  co 
o 
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-2 

which  decays  as  u  as  co  ->  ±°°;  this  spectral  decay  is  the  key 
issue  for  avoiding  a  logarithmic  singularity  in  e(t),  yH(t), 
y+(t),  and  y(t).  All  values  of  $  other  than  ±n/2  lead  to 
asymptotic  decay  of  Y(f)  in  (49)  according  to  -i  cos(4>)  co-1, 
which  leads  to  a  logarithmic  singularity  in  the  various  time 
functions  considered  here,  including  the  complex  envelope. 

Continuing  this  special  case  of  $  =  -  n/2  in  (67)  and  (68), 
we  find,  from  ( 48 ) , 

Z(f)  =  a'Tlco  '  z(t)  =  U(t>  (_i)  exP( -at )  for  t  f  •  (69) 

Also,  there  follows  from  (56),  (58),  and  (60),  respectively,  the 
error,  the  Hilbert  transform,  and  the  complex  envelope,  as 

e ( t )  =  -  £  Im(exp(-ct)  E^-ct))  ,  (70) 

yH(t)  =  -  U(t)  exp(-at)  cos(u>ot)  +  e(t)  ,  (71) 

y(t)  =  U(t)  (-i)  exp(-at)  +  i  exp(-iwQt)  e(t)  ,  (72) 

all  for  <{>  =  -  n/2. 

The  asymptotic  behavior  of  error  e(t)  at  infinity  is 
available  from  (3;  5.1.51]  as 


The  origin  behavior  is  available  from  [3;  5.1.11]: 


w 


e(t)  ~ 


o 


a 


+  co 


J. 

nt 


— z  as  t  ->  ±°° 


for  =  - 


o 
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e(t)  ~ 


1  - 


-  arctan(coQ/a) 


—  arctan(coo/a) 


as  t  ->  0- 


as  t  -»  0+ 


}  for  $  =  -  ~  .  (74) 


Observe  that  these  limits  in  (74)  at  t  =  ±0  are  both  finite. 

Also,  note  the  very  slow  decay  in  (73),  namely  1/t,  of  error  e(t) 
at  infinity. 

When  <|>  *  ±ji/2,  the  generalizations  to  (73)  and  (74)  are 
[3;  5.1.51  and  5.1.11] 


e(t)  ~ 


a  cos<j> 


to  sini 
o  T 


a2  + 


to. 


I 

nt 


— r  as  t  ->  ±® 


(75) 


and 


e(t)  ~  In | t |  as  t  ■>  0  . 


(76 


Now,  error  e(t)  becomes  infinite  at  the  origin  and  decays  only  as 
1/t  for  large  t.  (If  tan$  =  a/w0'  then  e(t)  =  0(t-2)  as  t  -*  ±«>; 
this  corresponds  to  Y(0)  =  0  in  (49).) 
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GENERAL  HILBERT  TRANSFORM  BEHAVIOR 

The  example  of  y(t)  in  (45)  (when  <}>  *  ±ji/2)  illustrates  the 
general  rule  that  if  a  time  function  has  a  discontinuity  of  value 
D  at  time  tQ,  then  its  Hilbert  transform  behaves  as  D/n  ln|t-tQ| 
as  t  •»  t  .  To  derive  this  result,  observe  that 

y(t)  ~  V  +  sgn( t-tQ )  as  t  UQ  ,  (77) 

when  y(t)  is  discontinuous  at  tQ.  Then,  for  t  near  t  ,  the 
Hilbert  transform  of  y(t)  is  dominated  by  the  components 


-  z 


...  1  f  du 

~  n  J  “ 


V  +  |d  sgn( t-tQ-u ) 


-b 

b 


+  n  I  ^  [V  +  5°  sgn(t-t0-U)j  , 


(78 


where  e  is  a  small  positive  quantity  and  the  principal  value 
nature  of  the  Hilbert  transform  integral  has  been  utilized.  The 
integrals  involving  constant  V  cancel;  also,  by  breaking  the 
integrals  in  (78)  down  into  regions  where  sgn  is  positive  versus 
negative,  and  watching  whether  t-t  is  positive  or  negative,  the 
terms  involving  ln(E)  cancel,  leaving  the  dominant  behavior 

yH(t)  ~  n  lnlt_t0l  as  t  *  to  ’  (79) 


(The  example  in  (66)  corresponds  to  a  discontinuity  D  =  cos(<j>)  at 
t  =  0,  as  may  be  seen  by  referring  to  (45).)  When  Hilbert 
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transform  VH(t)  has  this  logarithmic  singularity  (79),  then  so 
also  do  y+(t),  y(t),  and  e(t)  at  the  same  time  location.  Thus, 
the  complex  envelope  corresponding  to  a  discontinuous  y(t)  has  a 
logarithmic  singularity. 

An  alternative  representation  for  Hilbert  transform  yH(t)  in 
( 15 )  is  given  by 

yH(t)  =  |  df  exp( i2nft )  (-i)  sgn(f)  Y(f)  .  (80) 

If  Y ( f )  decays  to  zero  at  f  =  ±°°  and  if  Y  ( f )  is  continuous  for 
all  real  f,  then  an  integration  by  parts  on  (80)  yields  (due  to 
the  discontinuity  of  sgn(f))  the  asymptotic  decay 

yH(t)  ~  as  t  •*  ±«  .  (81) 

(Results  (73)  and  (75)  are  special  cases  of  (81),  when  applied  to 

example  (49).)  The  only  saving  feature  of  this  very  slow  decay 

for  large  t  in  (81)  is  that  Y(0)  may  be  small  relative  to  its 

maximum  for  f  *  0.  For  example  (49),  | Y ( f Q ) |  =  (2a)  *  for 

a  <<  wQ,  which  is  then  much  larger  than  Y(0)  =  -  sin(J>/wQ.  In 

this  narrowband  case,  the  slow  decay  of  (81)  will  not  be  overly 

significant  in  analytic  waveform  y+(t)  until  t  gets  rather  large. 

If  Y(0)  is  zero,  the  dominant  behavior  is  not  given  by  (81),  but 

2 

instead  is  replaced  by  a  1/t  dependence,  with  a  magnitude 
proportional  to  Y'(0). 
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GRAPHICAL  RESULTS 

We  now  take  the  example  in  (45)  with  parameter  values 
a  =  1  sec-*  and  f  =  1GU  Hz.  The  error  e(t)  in  (56)  is  plotted 
versus  time  t  in  figure  3  for  three  different  values  of  phase  <{> . 

A  time  sampling  increment  Afc  of  .02  msec  was  used  to  compute 
(56),  since  these  error  functions  are  very  sharp  in  t,  being 
concentrated  around  t  =  0  where  the  waveform  y(t)  has  its 
discontinuity.  The  period  of  the  carrier  frequency  is  1/f  =  10 

msec;  however,  the  error  functions  vary  significantly  in  time 
intervals  less  than  1  msec.  These  functions  approach  at 
t  =  0,  according  to  (66),  except  for  41  =  -n/2. 

The  corresponding  complex  envelope  is  given  by  (72);  its 
magnitude  is  plotted  in  figure  4  over  a  much  wider  time  interval. 
The  straight  line  just  to  the  right  of  the  origin  is  the  desired 
exponential  decay  a(t)  =  exp (-at),  which  dominates  the  error  e(t) 
in  this  region  of  time.  Eventually,  however,  for  larger  t  or 
negative  t,  the  error  e(t)  dominates,  with  its  much  slower  decay 
rate.  It  is  readily  verified  that  the  asymptotic  behavior 
predicted  by  (75)  is  in  control  and  very  accurate  near  both  edges 
of  figure  4. 

At  the  transition  between  the  two  components,  the  random 
vector  addition  in  (72)  leads  to  large  oscillations;  the  period 
of  the  carrier  is  1/f  =  1°  msec,  meaning  that  the  transition 

oscillations  in  figure  4  have  been  grossly  undersampled  with  the 
time  increment  A^  approximately  40  msec  that  was  used.  The  error 
curve  for  <j>  =  0  is  much  smaller  than  the  other  two  examples  over 
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most  of  its  range;  however,  the  magnitude  error  goes  sharply  to  <® 
at  t  =  0. 

From  (72)  and  figure  4,  it  is  seen  that  for  <j>  =  -n/2,  the 
phase  P(t)  of  complex  envelope  y.(t)  is  essentially  -n/2  for 
t  >  0,  until  we  reach  the  transition.  To  the  right  of  the 
transition,  the  phase  of  x(t)  exp(iu>ot)  is  essentially  n/2 
because  e(t)  >0  for  t  >  0,  for  this  example.  For  t  <  0,  the 
phase  of  ^(t)  exp(iu>ot)  is  -n/2  because  e(t)  <  0  for  t  <  0.  We 
will  numerically  confirm  these  claims  later  when  we  compute  the 
analytic  waveform  and  complex  envelope  by  means  of  FFTs. 


27 


TR  8827 


Figure  3.  Error  e(t)  for  Various  <j> 


Figure  4. 
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FILTERED  COMPLEX  ENVELOPE 

It  was  shown  in  (26)  that  the  spectrum  Y(f)  of  the  complex 
envelope  y(t)  °f  a  given  waveform  y(t)  with  complex  imposed 
modulation  z(t)  is  given  by  a  desired  term  Z ( f )  plus  an  undesired 
error  term,  namely. 


Y  (  f )  =  Z ( f )  +  i  E(f+f  )  .  (82) 

According  to  figures  1  and  2,  the  major  contribution  of  the  first 
term,  Z(f),  is  centered  around  f  =  0,  while  the  undesired  second 
term  in  (82)  is  centered  about  f  =  -f  .  This  suggests  the 
possibility  of  lowpass  filtering  complex  envelope  spectrum  Y(f) 
in  order  to  suppress  the  undesired  frequency  components.  Also, 
this  will  eliminate  or  suppress  the  undesired  logarithmic 
singularities  present  in  the  complex  envelope  y(t). 


LOWPASS  FILTER 

To  this  aim,  let 
and  cutoff  frequency 
filter  is  characteri 

H(f)  = 

The  filtered  complex 


H(f)  denote  a  lowpass 
f i ,  smaller  than  f  . 
ed  by 

cos2 |— |  for  | f |  < 

0  otherwise 
envelope  spectrum  is, 
G(f)  =  Y ( f )  H ( f  )  . 


filter  with  H(0)  =  1 
For  example,  the  Hann 

fl' 

-  •  (83) 

i 

in  general, 

(84) 
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The  importance  of  having  ^  <  f  is  that  filter  H(f)  will  then 
smoothly  cut  off  its  response  before  reaching  the  discontinuity 
at  f  =  -f  of  the  spectrum  Y(f)  of  the  complex  envelope  y(t);  see 
(27).  In  this  way,  we  can  avoid  the  slowly  decaying  behavior  of 
the  complex  envelope  y_{t)  for  large  t,  namely  1/t,  which 
inherently  accompanies  its  discontinuous  frequency  spectrum. 

This  will  prove  important  when  we  numerically  evaluate  the 
filtered  complex  envelope,  by  sampling  (84)  at  equispaced 
frequencies  and  performing  a  Fourier  transform  into  the  t  domain, 
necessarily  encountering  the  unavoidable  aliasing  in  time 
associated  with  such  a  technique. 

Since  the  complex  envelope  ^(t)  is  given  by  (25)  as  the  sum 
of  desired  component  z(t)  and  an  error  term,  the  filtered 
waveform  corresponding  to  spectrum  G(f)  in  (84)  is  given  by 

g(t)  =  y(t)  ©  h(t)  = 

=  z(t)  ©  h(t)  +  [i  e(t)  exp( -i2nfQt ) ]  ©  h(t)  =  (85) 

=  9d(t)  +  gu(t)  ,  (86) 

where  ©  denotes  convolution,  h(t)  is  the  impulse  response  of 
the  general  filter  H(f)  in  (84),  and  gd(t)  and  9u(t)  are, 
respectively,  the  desired  and  undesired  components  of  the 
filtered  complex  envelope  g(t).  We  should  choose  filter  H(f)  to 
be  real  and  even;  then  impulse  response  h(t)  is  also  real  and 
even . 

The  Hann  filter  example  in  (8?)  could  be  replaced  by  a  filter 
with  a  flatter  response  about  f  =  0  and  a  sharper  cutoff 
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behavior.  The  major  features  that  filter  H(f)  should  have  are  a 
fairly  flat  response  in  the  7(f)  frequency  range  near  f  =  0,  but 
cut  off  significantly  before  getting  into  the  major  frequency 
content  of  error  term  E(f+f  ),  which  is  centered  about  f  =  -fQ. 

If  the  given  waveform  y(t)  in  (1)  is  not  really  narrowrcnd,  there 
may  not  be  any  good  choice  of  cutoff  frequency  f ^ ;  that  is,  it 
may  be  necessary  to  sacrifice  some  of  the  higher  frequency 
content  of  z(t)  or  to  allow  some  of  the  error  e(t)  to  pass. 


EXAMPLE 


We  again  consider  the  example  given  in  (47)  and  (48),  along 
with  the  Hann  filter  in  (83).  In  order  to  evaluate  the  filtered 
complex  envelope  g(t)  in  (86),  we  define  an  auxiliary  function 

E ( z )  =  exp(z)  Ej(z)  ,  (87) 

where  E^(z)  is  the  exponential  integral  [3;  5.1.1].  Then,  when 
we  use  the  fact  that  (83)  can  be  expressed  as 

H  (  f  )  =  ^  exp(inf/f^)  +  ^  exp  (  -  in  f /f  ^  )  for  |f|  <  f^  ,  (88) 

we  encounter  the  following  two  integrals.  First,  we  need  the 
result 


£ 


-f 


exp  (  iwt+iorln/f  .  ) 
df 


a  +  i2w  +  iw 
o 


iz. i>" 

i  2  n 


exp(  -iujt.  )E(un)  -exp(  iu1t)E(vn) 


1 


(  89 
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where  co  =  2nf,  n  is  an  integer,  f^  <  f  ,  and  we  defined 

un  =  -(a  +  i2wQ  -  iwj)(t  +  ^n/fj)  , 

vn  =  -(a  +  i2w0  +  +  in/fi>  *  (90> 


To  derive  this  result,  we  let  x  =  -(a  +  i2wQ  +  ioo)(t  +  hn/f^)  in 
(89)  and  used  [3;  5.1.1],  along  with  the  important  fact  that 
ff  <  f  ,  which  guarantees  no  crossing  of  the  negative  real  avis 
of  the  resulting  contour  of  integration  in  the  complex  x-plane. 
Also,-  when  we  define 

-n  s  un(fo=0)  =  _(a  "  iwi)(fc  +  In/fi>  ' 

-n  H  vn(fo=0)  =  _(a  +  iwi)(t  +  \n/fl)  '  (91) 


then  for  fQ  =  0,  we  find  the  second  Integral  result  required, 
namely , 


exp(  icot  +  iw^-n/f . 
df  — 


a  +  ico 


izir 

i  2  n 


exp(-iw1t)E(un)-exp(iw1t)E(vn) 


1 


+  U(t  +  ^-n/fj)  exp(-at  -  ^-an/f^ 


(92) 


The  extra  term  in  the  second  line  of  (92)  is  due  to  a  crossing  of 
the  negative  real  axis  in  the  complex  x-plane  by  the  contour  of 
integration  when  we  make  the  substitution 
x  =  -  (a  +  iw)(t  +  ’jn/fj)  in  the  integral  of  (92). 

The  desired  component  of  the  filtered  complex  envelope  is 
given  by  the  first  term  of  (85)  and  (86),  in  the  alternative  form 
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gd(t)  =  J  df  exp(i2nft)  Z(f)  H(f)  = 


=  J  df  exp( i2nft )  §-Xg(^  cos2(s  Tt)  = 


(  93 


=  e 


i*(e“at[jU(t)  +  ±u(t  ♦  3^)  exp(^f-)  +  |u(t  -  ji-)  exp^f-) 
+  i^(exp{-i^1t)[E(!io)  -  \  Efuj)  -  i  B(u  _1i]  - 


-  exp(iw1t) 


E(v0)  -  ~  E(Yl)  -  \  E (v  _1) 


(94 


Here,  we  also  used  (48),  (88),  and  (92).  Since  the  factor 
multiplying  exp(i<j>)  in  (93)  has  conjugate  symmetry  in  frequency 
f,  the  time  function  multiplying  exp(i$)  in  (94)  is  purely  real 
for  all  time  t. 

The  undesired  spectral  component  in  (82)  is  given  by 


i  E(f+f  )  =  Z*(-f-2f  )  for  f  >  -f  , 
'  o  o  o 


(95) 


where  we  used  (22).  Therefore,  using  restriction  fj  <  fQ ,  the 
undesired  time  component  in  the  filtered  complex  envelope  in 
(86)  is  given,  upon  use  of  (89),  by 


gu(t: 


*  J  d£  exP<i2"ft>  a  cos2(i  It)  = 

-f,  ° 


-  e"PUnt)(elIP<-lMlt>[B(u0>  -  1  E(ul»  -  7  E(u-1> 


-  expjiWjt)  E( Vq )  -  3  E(Vj)  -  3  E(v_j)j)  .  (96 
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In  contrast  with  (94),  the  time  function  multiplying  phase  factor 
exp(-i<$>)  in  (96)  is  complex.  The  total  time  waveform  at  the 
filter  output,  g(t),  namely  the  filtered  complex  envelope,  is 
given  by  (94)  plus  (96),  and  depends  on  <{>.  In  fact,  since  the 
magnitude  of  total  output  g(t)  depends  on  we  will  look  at 
plots  of  the  magnitudes  of  components  |gd(t)|  and  |g  (t)|, 
neither  of  which  depend  on  <$> . 

For  comparison,  the  complex  envelope  itself  is  given  by  (26) 
in  the  form  y;(t)  =  z(t)  +  i  e(t)  exp( -i2nfQt ) .  Since  these  two 
(unfiltered)  components  depend  differently  on  phase  <j>,  we  shall 
also  consider  only  their  magnitudes  |z(t)|  and  |e(t)|  and  compare 
them  with  filtered  components  |gd(t)|  and  |gu(t)|,  respectively. 
In  particular,  from  (48),  the  desired  component  of  the  complex 
envelope  y(t)  for  the  example  at  hand  is 

z(t)  =  exp(i<)>  -  at)  U(t)  for  all  t  ,  (97) 

while  the  undesired  portion  is  given  by  (56)  and  (87)  as 

e  ( t )  =  -  ^Re{exp(i<j>  -  ct)  E^-ct)}  =  -  ^Re  j  exp(  i<j> )  E(-ct)}  (98) 

for  t  *  0,  where  c  =  a  -  iwQ.  The  magnitude  of  complex  waveform 
z(t)  is  independent  of  <|>,  but  the  magnitude  of  real  error 
waveform  e(t)  still  depends  on  <)>;  see  figure  3. 

The  magnitudes  of  z(t)  and  gd(t)  for  a  =  1  sec  *  and 
fj  =  40  Hz  are  displayed  in  figure  5  on  a  logarithmic  ordinate. 
The  filtered  complex  envelope  component,  gd(t),  drops  very 
quickly  to  the  left  of  t  =  0  and  is  indistinguishable  from  z(t) 
for  t  >  0;  compare  with  figure  4.  Thus,  the  passband  of  the  Hann 


34 


TR  8827 


filter  H ( f )  in  (83)  has  been  taken  wide  enough  to  pass  the 
majority  of  the  frequency  components  of  desired  function  Z(f)  in 
this  example.  The  darkened  portion  of  the  plot  just  to  the  left 
of  t  =  0  corresponds  to  a  weak  amplitude-modulated  40  Hz 
component,  which  is  the  cutoff  frequency  f^  of  filter  H(f). 

The  magnitudes  of  gu(t)  and  error  e(t)  are  displayed  in 
figure  6  for  the  additional  choice  of  parameters  f  =  100  Hz  and 
<{>  =  -n/2  rad.  The  peak  values  of  these  undesired  components  at 
t  =  0  differ  by  over  a  factor  of  10,  through  this  process  of 
filtering  the  complex  envelope.  At  the  same  time,  the  skirts  of 
filtered  version  gu(t)  are  down  by  several  orders  of  magnitude 
relative  to  e(t).  The  thick  plot  of  |gu(t)|  is  again  a  40  Hz 
component,  which  has  been  sampled  at  a  time  increment  =  .002 
sec. 

For  phase  <J»  =  0  instead,  original  waveform  y(t)  in  (45)  is 
discontinuous  at  t  =  0,  giving  rise  to  a  Hilbert  transform  which 
has  a  logarithmic  infinity  there;  see  (77),  (78),  and  (79). 
Therefore,  the  magnitude  of  error  e(t)  in  figure  7  has  an 
infinity  at  t  =  0,  whereas  the  filtered  quantity  9u(t)  is  finite 
there;  in  fact,  |g  (t)|  is  independent  of  <{> .  Although  e(t)  is 
significantly  reduced  in  value,  away  from  the  origin,  relative  to 
figure  6,  it  is  still  larger  than  the  filtered  quantity  gu(t). 
Since  the  energy  in  error  waveform  e(t)  is  independent  of  <j>  (see 
(62)),  smaller  skirts  in  e(t)  can  only  be  accompanied  by  a  larger 
peak;  in  fact,  this  latter  case  for  e(t)  in  figure  7  has  an 
infinite  (integrable)  peak  at  t  -  0.  By  contrast,  the  energy  in 
the  filtered  undesired  component  9u(t)  is,  from  (82)  and  (84), 
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J  df  | E ( f +f 0 ) | 2  I H ( f ) ! 2  ,  (99) 

which  can  be  considerably  less  than  the  error  energy,  when  filter 
H ( f )  significantly  rejects  the  displaced  error  spectrum  E(f+f  ). 

This  example  points  out  that  considerable  reduction  of  the 
undesired  error  term  in  the  complex  envelope  can  be  achieved 
through  the  use  of  lowpass  filtering  with  an  appropriate  cutoff 
frequency,  and  that  the  undesired  singularities  can  be  signif¬ 
icantly  suppressed.  Furthermore,  the  desired  component  of  the 
complex  envelope  can  be  essentially  retained.  These  conclusions 
follow  if  the  bandwidth  of  the  imposed  modulation,  z(t)  in  (1) 
and  (2),  is  small  relative  to  the  carrier  frequency  fQ. 
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TRAPEZOIDAL  APPROXIMATIONS  TO  ANALYTIC  WAVEFORM, 

COMPLEX  ENVELOPE,  AND  FILTERED  COMPLEX  ENVELOPE 

In  this  section,  we  address  methods  of  evaluating  the 
analytic  waveform  and  the  complex  envelope  by  means  of  FFTs . 

We  start  by  repeating  the  results  in  (6)  and  (8)  for  the 
analytic  waveform,  that  is, 

Y+(f)  =  2  U(f)  Y ( f )  ,  (100) 

oo 

y+(t)  =  J  df  exp(i2nft)  Y+(f)  =  J  df  exp(i2nft)  2  Y(f)  .  (101) 

0 

The  trapezoidal  approximation  to  (101)  is  obtained  by  sampling 
with  frequency  increment  A  to  get 

OO 

y+(t)  =  A  XZI  E  exp( i2nnAt )  2  Y(nA)  «  (102) 

n=0  n 

OO 

=  J  df  exp(i2nft)  2  Y(f)  A  &^(f)  = 

0 

=  y+(t)  ®  s1/A(t)  =  EZ  y+ft  -  g)  ,  (io3) 

where  sequence  eQ  =  h  and  en  =  1  for  n  >  1,  and  summations 
without  limits  are  from  -<*>  to  +°°. 

Notice  that  approximation  y+(t)  is  a  continuous  function  of 
time  t  and  has  period  1/A  in  t.  The  desired  term  in  (103)  is 
that  for  n  =  0,  namely  analytic  waveform  y+(t).  Because  y+(t) 
can  contain  a  slowly  decaying  Hilbert  transform  component,  the 
aliasing  at  separation  1/A  in  (103)  can  lead  to  severe  distortion 
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in  approximation  y+(t)  defined  in  (102). 

Since  y+(t)  has  period  1/A  in  t,  we  can  confine  its 
computation  to  any  interval  of  length  1/A.  In  particular,  if  we 
divide  this  interval  into  N  equally-spaced  points  (where  integer 
N  is  arbitrary),  we  can  compute,  from  (102), 

y+(jjf)  =  A  £11  en  exp(  i2nnk/N)  2  Y(nA)  (104) 

for  any  N  contiguous  values  of  k.  If  we  choose  the  range 
0  <  k  £  N-l,  and  if  we  collapse  the  infinite  sequence  in  the 
summand  of  (104)  according  to 


zn  =  2 A  En+jN  Y(nA  +  jNA)  for  0  <  n  <  N-l  , 


(105) 


then  (104)  can  be  written  precisely  as 


„  f  k)  Nzi. 

Y  + 1 NA  =  exp( i2nnk/N )  zn 

v  n=0 


( 106 


This  last  result  can  be  accomplished  by  means  of  an  N-point  FFT 
if  N  is  highly  composite.  This  is  a  very  efficient  method  of 
computing  the  aliased  version  of  the  analytic  waveform  as  defined 
by  ( 102  )  . 
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COMPLEX  ENVELOPE 

The  center  frequency  f  of  single-sided  spectrum  Y+(f)  in 
(100)  can  be  found  by  the  method  described  in  appendix  A.  Then 
the  complex  envelope  spectrum  and  waveform  are,  respectively, 

Y(f)  =  Y+(f+fc)  ,  (107) 

y(t)  =  J  df  exp(i2nft)  Y(f)  = 

=  |  df  exp( i2nf t )  Y+(f+fc)  «  exp(-i2nfct)  y+(t)  .  (108) 

The  approximation  to  complex  envelope  y^(t)  is  achieved  by 
relating  it  to  that  for  analytic  waveform  y+(t)  according  to 

y(t)  5  exp(-i2nfct)  y+(t)  =  (109) 

00 

=  exp(-i2nf  t)  A  >  '  e  exp(i2nnAt)  2  Y(nA)  ,  (110) 

c  n=0  n 

where  we  used  (108)  and  (102).  The  continuous  function 
exp(i2nfct)  £(t) ,  which  is  just  y+(t),  has  period  1/A  in  t,  which 
simplifies  its  calculation.  Using  (109),  (103),  and  (108),  there 
follows,  for  the  approximation  to  the  complex  envelope, 

Y(t)  =  exp(-i2nfct)  HZ  Y+(t  -  §)  =  HZ  xft  "  j)  exp( -i2nfcn/A)  . 

(Ill) 

The  desired  term  in  (111),  for  n  =  0,  is  complex  envelope  ^(t). 
The  n-th  term  has  a  time  delay  (aliasing)  of  n/A  and  a  phase 
shift  of  n2nfc/A  radians,  which  is  arbitrary  because  frequency 
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sampling  increment  A  in  (102)  is  unrelated  to  center  frequency 
fc  of  Y+(f)  in  (100). 

Sample  values  of  complex  envelope  approximation  y(t)  can  be 
obtained  from  (110)  as 


exp 


-i2nf 


_k 
C  NA 


en  exp(i2nnk/N) 


2  Y(nA) 


(112) 


Again,  the  infinite  sum  in  (112)  can  be  collapsed  and  realized 
as  an  N-point  FFT ;  see  (104)  -  (106).  The  phase  factor 
Pk  =  exp( -i2itfck/(NA)  )  can  be  computed  via  recurrence 
pk  =  Pk_i  exp(-i2nfc/(NA) ) . 


FILTERED  COMPLEX  ENVELOPE 

The  spectrum  of  the  filtered  complex  envelope  is  given  by 
(84)  as  G(f)  =  Y(f)  H(f).  The  filtered  complex  envelope  waveform 
is 

g ( t )  =  |  df  exp( i2nft )  G(f)  =  y(t)  ©  h(t)  (113) 

and  has  low  sidelobes  and  rapid  decay  in  t,  when  filter  H(f)  is 
chosen  appropriately. 

The  approximation  to  g(t)  adopted  here  will  be  generalized 
slightly  in  order  to  allow  for  frequency-shifted  sampling. 
Specifically,  we  define 
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gQ(t)  m  J  df  exp( i2nft )  G(f)  A  &A(f  -  a)  =  (114) 

=  A  )  exp(i2n[nA  +  a]t)  G(nA  +  a)  .  (115) 

n 

The  function  exp(-i2nat)  gQ(t)  has  period  1/A  in  t,  which  allows 
us  to  confine  its  calculation  to  any  convenient  period. 

The  behavior  of  approximation  gQ(t)  in  (114)  follows  as 

gQ(t)  =  g(t)  e  [exp(i2nat )  &1/A(t)]  = 

=  g(t)  ©  )  '  exp(i2nan/A)  &(t  -  = 

n  1 

=  )  '  exp( i2nan/A)  gft  -  £]  .  (116) 

This  is  the  aliased  version  of  the  filtered  complex  envelope. 

The  desired  term,  for  n  =  0,  is  the  filtered  complex  envelope 
g(t),  independent  of  the  choice  of  frequency  shift  a.  Shift  a 
is  arbitrary  and  could  be  taken  as  -f  if  desired. 

Samples  of  ga(t)  are  available  from  (115)  according  to 

(nS)  =  A  exp(i2na  2ZI  exp(i2nnk/N)  G(nA  +  a)  ,  (117) 

which  we  can  limit  to  0  $  k  1  N-l  due  to  the  periodicity  of 
ga(t).  Again,  the  infinite  sum  on  n  can  be  converted  to  an 
N-point  FFT  without  error,  by  collapsing  into  the  finite  sequence 

z  =  A  >  G(nA  +  a  +  jNA)  for  0  <  n  <  N-l  .  (118) 

n  j 
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The  remaining  phasor  exp( i2nak/( NA ) )  in  (117)  can  be  quickly 
obtained  via  recursion  on  k. 


GRAPHICAL  RESULTS 

The  same  fundamental  example  introduced  in  (45)  will  be 
used  here,  again  with  a  =  1  sec-*  and  f  =  100  Hz.  For  phase 
<)>  =  -n/2,  FFT  size  N  =  1024,  and  a  frequency  increment  of 
A  =  1/80  Hz,  the  magnitude  of  y/t),  namely  A(t),  is  displayed  in 
figure  8  over  the  1/A  =  80  sec  period  centered  at  t  =  0 .  This 
selection  of  the  time  period  has  been  purposely  made  the  same  as 
that  used  in  figure  4,  for  easy  comparison  of  results.  The  major 
difference  between  the  <)>  =  -n/2  result  in  figure  4  and  figure  8 
is  that  the  aliasing  in  the  latter  case  causes  the  curve  to  have 
a  jagged  behavior  and  to  droop  in  the  neighborhood  of  t  =  ±40 
sec.  However,  other  examples  could  well  have  the  aliasing 
increase  near  the  edges  of  the  period.  A  total  of  88,000  samples 
of  Y ( f )  at  frequency  increment  A  were  taken  in  computation  of 
(104);  the  collapsing  in  (106)  resulted  in  storage  of  only 
N  =  1024  complex  numbers  and  the  ability  to  use  a  single 
relatively  small  N-point  FFT.  A  program  for  the  evaluation  of 
the  complex  envelope  by  means  of  an  FFT  with  collapsing  is 
furnished  in  appendix  B;  the  FFT  uses  a  zero-subscripted  array  in 
direct  agreement  with  the  mathematical  definition  of  the  FFT. 

The  corresponding  phase,  P(t)  =  arg{y(t)j,  of  the  aliased 
complex  envelope  is  given  in  figure  9.  The  phase  is 
approximately  -n/2  for  0  i  t  1  10  sec,  as  expected,  since  in 
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this  limited  time  interval,  the  error  is  not  the  dominant  term. 
However,  over  the  rest  of  the  period,  the  error  term  does 
dominate  and  it  has  an  exp(-iwQt)  behavior,  where  f  =  100; 
see  (25)  and  figure  3.  Thus,  the  time  sampling  increment 
=  1/ ( NA )  =  .078  sec  is  grossly  inadequate  to  track  this 
high-frequency  term,  and  we  get  virtually  random  samples  of  the 
phase  of  the  complex  exponential  exp( -iwQt ) . 

To  confirm  the  phase  behavior  outside  the  (0,10)  sec 
interval,  we  have  plotted  the  phase- t>f  y(t)  exp(iwQt)  =  y+(t)  in 
figure  10  as  found  by  the  FFT  procedure  above.  To  the  right  of 
t  =  10,  the  phase  is  approximately  n/2,  in  agreement  with  the 
fact  that  e(t)  is  real  and  positive  for  t  >  0;  see  figure  3  and 
(25).  For  time  t  <  0,  the  phase  is  -n/2  because  e(t)  <  0  for 
t  <  0.  The  oscillatory  behavior  at  both  edges  of  the  period, 
namely,  for  30  <  |t|  <  40,  is  due  to  aliasing  from  adjacent  lobes 
indicated  by  (103)  and  (111). 

When  4>  is  changed  to  0  and  everything  else  is  kept  unchanged, 
the  result  for  the  magnitude  of  complex  envelope  aliased  version 
y(t)  is  plotted  in  figure  11.  Comparison  with  the  exact  results 
in  figure  4  reveals  a  very  dramatic  increase  in  aliasing,  in 
fact,  by  two  orders  of  magnitude.  The  reason  for  this 
considerable  increase  can  be  seen  from  figure  3  and  (75);  namely, 
the  error  e(t)  is  unipolar  for  $  =  0  and  it  decays  very  slowly. 
Whereas  for  figure  8,  the  alternating  character  of  the  over¬ 
lapping  aliased  error  lobes  led  to  a  cancellation  near  t  =  ±40 
sec,  the  opposite  situation  occurred  in  figure  11,  leading  to  a 
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considerable  build-up  of  the  aliasing  effect. 

The  corresponding  phase  of  ;y(t),  P(t),  is  plotted  in  figure 
12.  Its  value  is  zero  in  the  region  0  d  t  <  10,  as  expected, 
since  the  desired  term,  exp(-at),  dominates  here.  Outside  this 
region,  the  situation  is  the  same  as  explained  above  with  respect 
to  figure  9.  We  have  not  plotted  the  counterpart  to  figure  10 
because  no  one  error  lobe  dominates  anywhere  on  the  time  scale; 
the  result  is  a  phase  plot  that  looks  random  over  the  entire 
period  of  (-40,40)  sec. 

When  the  complex  envelope  spectrum  is  filtered  according  to 
the  Hann  filter  in  (83)  -  (86),  the  results  for  the  sampled 
filtered  complex  envelope  waveform,  obtained  by  means  of  the 
collapsed  FFT  in  (117)  and  (118)  with  a  =  u,  are  given  in  figures 
13  and  14.  There  were  6400  frequency  samples  taken  uf  G(f)  with 
increment  A  =  1/80  Hz  and  an  FFT  size  of  N  =  1024  was  utilized; 
see  appendix  B.  A  comparison  of  the  magnitudes  in  figures  13  and 
5  reveals  virtually  identical  results;  namely,  the  error  and  its 
inherent  accompanying  aliasing,  that  was  present  in  figure  8,  is 
absent  from  figure  13. 

The  corresponding  phase  plot  of  the  FFT  output  is  displayed 
in  figure  14.  in  the  region  0  1  t  1  24  sec,  where  the  desired 
exp(-at)  term  dc  mates,  the  FFT  output  phase  is  equal  to  the 
value  of  <(>  =  -n/2  for  this  example.  When  this  example  was  rerun 
for  <f>  =  0 ,  similar  high  quality  results  were  obtained,  except 
that  the  FFT  output  phase  was  zero.  The  benefits  of  filtering 
the  complex  envelope  spectrum  are  well  illustrated  by  the  results 
of  figures  13  and  14. 
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Figure  8.  Magnitude  of  Complex  Envelope  via  FFT ,  $  =  - rt/2 


Figure  9.  Phase  of  Complex  Envelope  via  FFT,  y=-n/2 
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Figure  10.  Phase  of  Analytic  Waveform  via  FFT ,  <p  =  -n/2 
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Figure  11.  Magnitude  of  Complex  Envelope  via  FFT,  <J>  =  0 


Figure  12.  Phase  of  Complex  Envelope  via  FFT ,  <J>  =  0 
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ALIASING  PROPERTIES  OF  COSINE  AND  SINE  TRANSFORMS 

If  a  time  function  is  causal,  it  can  be  obtained  from  its 
Fourier  transform  either  by  a  cosine  or  a  sine  transform. 

However,  when  these  integral  transforms  are  approximated,  by 
means  of  sampling  the  frequency  function  and  using  some 
integration  rule  like  trapezoidal,  the  "alias-free"  interval  in 
the  time  domain  is  approximately  halved,  as  shown  below.  This 
does  not  necessarily  mean  that  these  transform  alternatives 
should  be  discarded,  because  a  more  rapidly  decaying  integrand 
can  be  useful,  but  it  does  point  out  a  cautionary  feature  in 
their  use  and  the  need  to  consider  the  tradeoff  between  aliasing 
and  truncation  error. 

GENERAL  TIME  FUNCTION 

In  general,  complex  time  function  y(t)  is  obtained  from  its 
Fourier  transform  Y(f)  according  to 

y(t)  =  J  df  exp(i2nft)  Y(f)  =  (119) 

=  J  df  cos(2nft)  Y ( f )  +  i  J  df  sin(2nft)  Y(f)  =  (120) 

■  Ye(t)  +  y0(t)  for  a11  t  /  (i2i) 

where  complex  functions  ye(t)  and  y o ( t )  are  the  even  and  odd 
parts  of  function  y(t),  respectively. 
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CAUSAL  COMPLEX  TIME  FUNCTION 

Now  suppose  that  y(t)  is  causal,  but  possibly  complex;  then 

y(t)  =  0  for  t  <  0  .  (122) 

Then,  letting  t  =  -a,  a  >  0,  we  have,  from  (122)  and  (121), 

0  =  y ( -a )  =  ye(~a)  +  yQ(-a)  =  Ye(a)  -  yQ(a)  for  a  >  0  .  (123) 

That  is, 

yo(a)  =  ye(a)  for  a  >  0  .  (124) 

Therefore,  from  (121)  and  (120),  we  have  two  alternatives  for  a 
causal  complex  time  function  y(t): 

y(t)  =  2  J  df  cos(2rtft)  Y(f)  for  t  >  0  ,  (125) 

and 

y(t)  =  i2  J  df  sin(2nft)  Y(f)  for  t  >  0  .  (126) 

We  need  complex  function  Y(f)  for  negative  as  well  as  positive 
frequency  arguments  f,  in  order  to  determine  causal  complex 
function  y(t),  but  we  can  utilize  either  a  cosine  or  a  sine 
transform . 
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NONCAUSAL  REAL  TIME  FUNCTION 

Now  suppose  instead  that  y(t)  is  real,  but  noncausal.  Then, 
since  spectrum  Y(-f)  =  Y*(f),  we  can  express  (119)  as 

00 

y(t)  =  2  Re  J  df  exp(i2nft)  Y(f)  =  (127) 

0 

00  CO 

=  2J  df  cos(2nft)  Yr(f)  -  2J  df  sin(2nft)  Y^(f)  for  all  t.  (128) 
0  0 

The  first  term  in  (128)  is  even  part  ye(t),  while  the  second  term 
in  (128)  is  odd  part  yQ(t);  see  (121).  In  this  case  of  a  real 
time  function  y(t),  we  need  complex  function  Y(f)  only  for  f  >  0. 

CAUSAL  REAL  TIME  FUNCTION 

Now  let  y(t)  be  both  causal  and  real.  Then  using  property 
Y(-f)  -  Y*(f)  in  (125)  and  (126),  we  obtain 

00 

4  J  df  cos(2nft)  Yr(f)  for  t  >  0  ,  (129) 

0 

00 

4  J  df  sin(2nft)  Y^(f)  for  t  >  0  .  (130) 

0 

Here,  we  need  either  Yr(f)  or  Y^(f),  and  then  only  for  positive 
frequency  arguments  f.  Also,  a  cosine  or  a  sine  transform  will 
suffice  for  determination  of  y(t). 


y(t)  = 


and 


y(t)  =  - 
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ALIASING  PROPERTIES 

The  above  relations  have  all  assumed  that  spectrum  Y(f)  is 
available  for  all  continuous  f.  Now  we  will  address  the  effects 
of  only  having  samples  of  Y(f)  available  at  frequency  increment 
A.  We  begin  with  the  trapezoidal  approximation  to  (119): 

y^(t)  =  A  Y1  exp(i2nnAt)  Y(nA)  for  all  t  .  (131) 

n 

The  approximation  y ^ ( t )  is  periodic  in  t  with  period  1/A.  It 


can  be  expressed  exactly  as 

y^(t)  =  J  df  exp(i2nft)  Y(f) 

A  &A(f)  - 

(132) 

=  y(t)  ®  s1/4(t)  =  E  y(t  -  §) 

n 

for  all  t  . 

(133) 

That  is,  approximation  y^(t)  is  an  aliased  version  of  desired 
waveform  y(t),  with  displacements  1/A  in  time.  This  result  holds 
for  any  complex  waveform  y(t)  and  has  been  used  repeatedly  in  the 
analyses  above. 

The  second  approximation  of  interest  is  obtained  from  the 
cosine  transform  in  (125),  which  applies  for  causal  complex  y(t) 
in  the  form 

; y 2  (t)  =  2A  ]F  cos(2nnAt)  Y ( nA )  for  all  t  .  (134) 

n 

y2c(t)  also  has  period  1/A  in  t  and  can  be  developed  as  follows: 
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y9  (t)  =  A  Y1  [exp(  i2ixnAt )  +  exp( -i2nnAt ) ]  Y(nA)  = 
zc  n 

-  J  df  [exp(i2nft)  +  exp(-i2nft)]  Y(f)  A  &^(f)  = 

=  [y(t)  +  y ( — t ) ]  ©  &1/A(t)  =  2  ye(t)  ©  &1/A(t)  = 

=  E  [y(t  -  |)  +  y(f  -  t)j  for  a11  t  .  (i35) 

That  is,  sampling  of  the  cosine  transform  in  (125)  results  in 
aliasing  of  y(t)  plus  its  mirror  image  y(-t),  even  when  y(t)  is 
causal.  This  will  restrict  useful  results  in  y2C(t)  to  a  region 
approximately  half  as  large  as  that  given  by  (131)  and  (133), 
where  the  sampled  exponential  transform  was  used.  Even  when  we 
restrict  calculation  of  approximation  y2C(t)  to  the  region 
(0,1/A),  we  are  contaminated  by  the  mirror  image  lobe  y(l/A  -  t) 
and  by  the  usual  lobe  y(t  +  1/A)  extending  from  t  =  -1/A  into  the 
desired  region. 

A  similar  situation  exists  for  using  a  sampled  version  of  the 
sine  transform  for  causal  complex  y(t)  in  (126);  namely,  consider 
the  approximation 

y9  (t)  s  i2A  Y1  sin(2nnAt)  Y(nA)  for  all  t  .  (136) 

n 

Then 


55 


TR  8827 


y2S(t)  =  A  E  [ exp( i2nnAt )  -  exp( -i2nnAt ) ]  Y(nA)  = 
n 

=  J  df  [exp'i2nft)  -  exp( -i2nft ) ]  Y(f)  A  &^(f)  = 


=  [y(t)  -  y ( -t) ]  ©  &1/A(t)  =  2  yo(t)  ©  &1/A(t)  = 


=  E[y(t  -  g)  -  y  (§  -  t)  ]  for  all  t  . 


(137) 


Here,  for  the  approximate  sine  transform,  twice  the  odd  part  of 
causal  complex  y(t)  is  aliased  with  separations  1/A  in  time, 
thereby  again  leading  to  a  clear  region  only  about  half  that 
attainable  from  (131)  and  (133).  We  will  return  to  these 
apparently  undesirable  transform  properties  below  and  find  them 
useful  when  we  consider  a  causal  real  time  function. 

The  next  approximation  is  for  the  noncausal  real  waveform 
result  in  (127);  namely,  letting  eq  =  h  and  en  =  1  for  nil,  we 
have  trapezoidal  approximation 


y ^ ( t )  =  2  Re  A  )  '  exp(i2nnAt)  Y(nA)  for  all  t  .  (138) 

n=0 


Then 


y^(t)  =  2  Re  J  df  exp(i2nft)  Y(f)  A  &A(f)  = 

0 

=  J  df  exp(i2nft)  Y(f)  A  &A(f)  =  E  y(fc  “  for  a11  t  •  (139) 
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just  as  in  (133).  Thus,  the  combination  of  the  cosine  and  sine 
transforms  in  (128)  does  not  additionally  damage  the  aliasing 
behavior  associated  with  sampling.  In  practice,  we  would  use  the 
real  part  of  the  exponential  transform  as  given  by  (127).  The 
same  result,  (139),  follows  when  the  cosine  and  sine  transforms 
in  (128)  are  individually  directly  approximated  by  the 
trapezoidal  rule  and  the  results  added  together. 

The  two  final  approximations  of  interest  come  from  sampling 
the  results  for  causal  real  y(t)  in  (129)  and  (130);  from  (129), 
define  approximation 

00 

y.  (t)  2  4A  )  '  e  cos(2nnAt)  Y  (nA)  for  all  t  ,  (140) 

qc  n=0  n  r 

which  has  period  1/A  in  t.  Now  we  develop  (140)  as 

y4c(t)  =  2A  51  cos( 2nnAt)  Yr(nA)  = 
n 

=  2  J  df  cos(2nft)  Yr(f)  A  &^(f)  = 

=  2  J  df  exp(i2nft)  Yr(f)  A  6^(f)  = 

=  J  df  exp( i2nft )  [Y(f)  +  Y*(f)]  A  6A ( f )  = 

=  y ( t )  ®  &1/A(t)  +  y*(-t)  ©  &1/A(t)  = 

=  H  y(t  *  f )  +  £  y(f  -  t)  =  2  ye(t)  ©  &1/A(t)  for  all  t,  (141) 
n  '  n  ' 
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where  we  used  the  real  character  of  y(t). 

This  end  result  is  identical  to  (135);  however,  approximation 
y4c(t)  in  (140)  uses  only  the  real  part  Yr(f)  of  the  complex 
function  Y(f),  whereas  Y2C(t)  i-n  (134)  requires  the  complete 
complex  function  Y(f)  for  a  causal  complex  y(t).  Since  it  is 
possible  to  have  complex  functions  Y(f)  which  have  rapidly 
decaying  real  parts  and  slower  decaying  imaginary  parts,  (140) 
affords  the  possibility  of  getting  a  smaller  truncation  error 
than  (134),  when  y(t)  is  causal  real  and  when  both  sums  are 
carried  out  to  the  same  frequency  limit,  because  both  sums  must 
be  terminated  in  practice.  Whether  the  reduction  in  the  usable 
"alias-free"  region,  dictated  by  (141),  can  be  traded  off  against 
a  smaller  truncation  error  associated  with  use  of  only  the  real 
part  Yr(f)  in  (140),  depends  on  the  particular  example  under 
investigation.  In  any  event,  (140)  affords  an  alternative  to 
consider  for  causal  real  y(t). 

The  final  approximation  comes  about  by  sampling  (130): 

CO 

y.  (t)  =  -4A  )  ‘  sin(2nnAt)  Y.(nA)  for  all  t  ,  (142) 

qs  n=l  1 

which  has  period  1/A  in  t.  In  the  usual  fashion,  we  find 
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y4s(t)  =  _2A  ^  sin( 2nnAt )  Yi(nA)  = 

=  -2  J  df  sin(2nft)  Y^(f)  A  &A(f)  = 

=  i2  J  df  exp(i2nft)  Y^(f)  A  &A(f)  = 

-  J  df  exp{ i2nft )  [Y(f)  -  Y* ( f ) ]  A  &A(f)  = 

-  y(t)  ©  *1/A( t)  -  y*(-t)  ©  &1/A( t)  =  2  yo(t)  ©  &1/A(t)  = 

=  £  yft  -  5)  -  £  y(|  -  t)  for  all  t  .  (143) 

Here,  we  used  the  real  character  of  y(t). 

The  end  result  in  (143)  is  identical  to  (137);  however, 
y4s(t)  in  (142)  only  requires  knowledge  of  the  imaginary  part 
Y^(f)  of  complex  function  Y(f),  whereas  y2S(t)  (136)  requires 
the  complete  complex  function  Y(f)  for  a  causal  complex  y(t). 

This  is  due  to  the  fact  that  (129)  and  (130)  apply  only  to  causal 
real  y(t),  whereas  (125)  and  (126)  apply  to  causal  complex  y(t). 
Since  there  exist  complex  functions  Y(f)  which  have  more  rapidly 
decaying  imaginary  parts  than  real  parts,  the  opportunity  arises 
to  reduce  the  truncation  error  by  employing  (142)  instead  of 
(136),  when  y(t)  is  causal  and  real.  The  comments  in  the  sequel 
to  (141),  regarding  the  trade-off  between  truncation  error  and  a 
reduced  alias-free  region,  are  again  applicable. 
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This  procedure,  of  using  only  the  imaginary  part  of  a  Fourier 
transform  because  it  decays  faster  than  the  real  part,  was 
utilized  to  advantage  in  [4;  pages  4-6]  and  was  based  upon  an 
earlier  result  in  [5;  (15)].  The  very  rapid  decay  of  the 
imaginary  part  far  outweighed  the  aliasing;  see  [4;  page  6]. 


EVALUATION  BY  MEANS  OF  FFTs 

If  periodic  function  y^(t)  in  (131)  is  evaluated  at  the 
equally  'paced  time  points  k/(NA)  for  k=0  to  N- 1 ,  which  suffice 
to  cover  one  period,  we  obtain 

yi(nI)  =  A  £  exp( i2nnk/N)  Y(nA)  =  (144) 

N-l 

_  ^  ^  *  exp(i2nnk/N)  z  ,  (145) 

n=0 

where  lzn)/  0  S  n  ^  N-l,  is  the  collapsed  version  of  sequence 
(Y(nA)},  -®  <  n  <  °°.  No  approximations  are  involved  in  this 
collapsing  procedure  from  (144)  to  (145).  Relation  (145)  can  be 
accomplished  by  means  of  an  N-point  FFT  if  N  is  highly  composite. 

In  a  similar  fashion,  (140)  yields  samples  of  the  cosine 
transform  as 
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y4c (nS)  =  4A  C  en  cos(2nnk/N)  Yf(nA)  = 

GO 

=  4A  Re  )  exp( i2nnk/N )  En  Yr(nA)  =  (146) 

n=0 

N-l 

=  4A  Re  )  '  exp( i2nnk/N)  z  ,  (147) 

n=0  n 

where  { zn } ,  0  i  n  £  N-l,  is  the  collapsed  version  of  sequence 
{ en  Yr(nA) } ,  0  £  n  <  ». 

Since  (147)  will  likely  be  realized  as  the  real  part  of  an 
FFT  output,  the  question  arises  as  to  the  interpretation  and 
utility  of  the  total  complex  FFT  output  in  (147).  To  this  aim, 
we  rewrite  y4c(t)  in  (146)  (in  its  continuous  time  version)  as 

<X> 

y4c(t)  =  Re  4  |  df  exp{ i2nft )  Yr(f)  A  &A( f )  = 

0 

=  Re { z x ( t )  ®  &1/A(t) )  ,  (148) 

where  we  define,  for  all  t,  Fourier  transform 

00 

z^t)  =  4  |  df  exp(i2nft)  Yr(f)  = 

0 

00 

=  I  df  exp( i2nft )  [2Y(f)  +  2Y*(f)]  = 

0 

=  Y+(t)  +  y*( -t)  =  [y(t)  +  y*(-t)]  exp(i2nfct)  = 

=  y ( t )  +  y ( -t)  +  i[yH(t)  -  yH(-t)]  .  (149) 
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That  is,  Y4C(t)  is  the  real  part  of  the  aliased  version  of  z^(t), 
which  itself  is  composed  of  the  analytic  waveform  y+(t)  and  its 
mirror  image.  Thus,  not  only  is  z^t)  aliased  according  to 

(148) ,  but  in  addition,  z^(t)  contains  terms  which  will  further 
overlap  and  thereby  confuse  the  values  of  y4c(t)  in  the 
fundamental  range  (0,1/A).  (Of  course,  the  real  part  of  z^(t)  in 

(149)  for  t  >  0  is,  as  expected,  just  y(t)  for  this  causal  real 
case .  ) 

Finally,  sampling  the  sine  transform  of  Y^(f)  in  (142)  yields 

,  °° 

y4s  NA  =  -  4A  2ZI  sin( 2nnk/N )  Yi(nA)  = 
v  '  n=l 

00 

=  -  4A  Im  )  '  exp(i2nnk/N)  Y.(nA)  =  (150) 

- r*  1 

n=l 

N-l 

=  -  4A  Im  )  '  exp{ i2nnk/N )  z  ,  (151) 

n=0  n 

where  (znl/  0  ^  n  *  N-l,  is  the  collapsed  version  of  sequence 
(Y^(nA),  1  <  n  <  “.  Relation  (151)  can  be  realized  as  an  N-point 
FFT  of  which  only  the  imaginary  part  is  kept  for  0  1  k  £  N-l. 

As  above,  the  interpretation  of  the  complete  complex  output 
of  the  FFT  in  (151)  is  furnished  by  returning  to  the  continuous 
version  of  the  sampled  y 4 s ( t )  in  (150).  We  express  it  as 

00 

y4s(t)  =  ~  Im  4  |  df  exp(  i2nf  t )  Y^f)  A  6^(f)  = 

0 

=  lm{z2(t)  ©  61/A(t) |  ,  (152) 
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where  we  define,  for  all  t,  Fourier  transform 

00 

z2(t)  =  -  4  J  df  exp(i2nft)  Y^(f)  = 

0 

00 

=  i  I  df  exp( i2nft )  (2Y(f)  -  2 Y* ( f ) ]  = 

0 

=  i[y+(t)  -  y* ( -t) ]  =  i[y( t)  -  **(-t)]  exp(i2nfct)  =  (153) 

=  i[y(t)  -  y(-t) ]  -  yH(t)  -  yH(-t)  .  (154) 

Again,  the  aliasing  of  z2(t)  in  (152)  and  the  mirror  image  of  the 
analytic  waveform  and  complex  envelope  in  (153)  will  serve  to 
confuse  the  usefulness  of  z2(t).  The  imaginary  part  of  z2(t)  in 
(154)  for  t  >  0  is  just  y(t),  as  expected,  for  this  causal  real 
waveform. 
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DISPLACED  SAMPLING 


If  displaced  samples  of  a  waveform  are  desired,  such  as  at 
time  locations  (k+0)/(NA)  in  (145),  where  0  <  £  <  1,  we  can 
obtain  them  via  an  N-point  FFT  as  follows:  from  (131), 


Ylp^f  =  A  exp(i2nnk/N)  exp(i2nn0/N)  Y(nA)  =  (155) 


N-l 

=  A  )  '  exp(i2nnk/N)  z  for  0  1  k  £  N-l  , 

n=0  n 


(156) 


where  { zn 1 ,  0  <  n  £  N-l,  is  the  collapsed  version  of  sequence 
{exp(  i2nn(5/N)  Y ( nA )  ) ,  -«>  <  n  <  ®.  That  is,  we  have  to  load  up 
the  arrays  containing  (zn)  with  phase-shifted  versions  of  the 
original  sequence  {Y(nA)j  and  then  perform  the  N-point  FFT. 
Calculation  of  phasor  pn  =  exp(i2nn0/N)  in  (155)  can  take 
advantage  of  recursion  pn  =  Pn-1  exp(i2nfS/N). 
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SUMMARY 

The  advantages  of  filtering  the  complex  envelope  spectrum  by 
means  of  a  suitable  lowpass  filter  are  significant  in  some 
instances.  The  singular  behavior  of  the  complex  envelope 
waveform  is  eliminated  by  utilizing  a  filter  which  cuts  off  at 
finite  frequencies,  while  the  slow  decay  in  the  time  domain  of 
the  complex  envelope  is  circumvented  by  using  a  filter  with  a 
smoothly  tapered  cutoff  that  prevents  any  discontinuities  in  the 
complex  envelope  spectrum  from  contributing. 

The  use  of  an  FFT  to  evaluate  the  filtered  complex  envelope 
is  then  an  attractive  efficient  approach  because  the  inherent 
time  aliasing  associated  with  frequency  sampling  has  been  greatly 
suppressed.  Also,  the  very  rapidly  varying  singular  components 
of  the  complex  envelope  have  been  eliminated,  allowing  for  a 
lower  time-sampling  rate,  that  is,  smaller  FFT  sizes. 

When  two  waveforms,  each  with  its  own  imposed  amplitude-  and 
phase-modulations,  are  convolved,  such  as  encountered  in  the 
narrowband  excitation  of  a  passband  filter,  the  output  complex 
envelope  is  given  exactly  by  the  convolution  of  the  individual 
complex  envelopes.  Although  the  convolution  of  the  two  (complex) 
imposed  modulations  is  often  a  good  approximation  to  the  output 
complex  envelope,  it  has  an  error  term.  This  analysis  is 
presented  in  appendix  C. 
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APPENDIX  A.  DETERMINATION  OF  CENTER  FREQUENCY 


Suppose  we  are  given  spectrum  Y(f)  of  (narrowband)  real 
waveform  y(t),  but  the  center  frequency  of  Y(f)  is  not  obvious  or 
is  unknown.  The  analytic  waveform  is  still  uniquely  given  by 

00 

y+(t)  =  J  df  exp(i2nft)  Y+(f)  =  2  J  df  exp(i2nft)  Y(f)  .  (A-l) 

0 

Make  a  guess  at  initial  frequency  f^  near  the  center  of  Y+(f). 
Then  compute  the  initial  down-shifted  waveform 


y^t) 


=  exp(-i2nf^t)  y+(t)  =  2  J  df  exp(i2nft)  Y(f+f 

-f . 


( A-2 


Compute  initial  phase  P^(t)  =  arg{y^(t))  and  then  unwrap  P^(t). 
Select  time  t  in  the  interval  T  of  interest  and  fit  a  straight 
line  a  +  pt  to  the  unwrapped  phase  P^(t)  over  T.  Compute 
frequency 


f 

c 


f  +  — £■  • 
l  2n  ' 


( A—  3 ) 


this  is  the  center  frequency  of  y+(t)  for  t  e  T.  Another 
selection  of  a  different  time  interval  could  lead  to  a  somewhat 
different  center  frequency;  there  is  no  unique  center  frequency 
of  an  arbitrarily  given  spectrum  Y(f). 

The  complex  envelope  is  then 


Y(t)  =  exp(-i2nfct)  y+(t)  .  (A-4) 

The  "physical"  envelope  or  extracted  amplitude  modulation  is 
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A(t)  =  I Y ( t ) |  =  |y+(t) !  =  |y.(t)|  ,  (A-5) 

which  is  independent  of  the  choice  of  f^  or  fc>  The  extracted 
phase  of  complex  envelope  ^{t)  is 

P(t)  «  arg  (y(t)  }  =  arg(y+(t)  exp ( -i2nf Qt ) }  =  Pi(t)  -  pt  ,  (A-6) 

where  we  used  (A-4)  and  (A-2).  Functions  y^(t)  and  P^(t)  have 
already  been  computed  and  can  be  used  to  evaluate  the  envelope 
A(t)  and  phase  P(t).  The  real  waveform  is 

y(t)  =  Re{y+(t) ]  =  Re{y(t)  exp(i2nfct)}  =  A(t)  cos[2nfct  +  P(t)], 

( A-7  ) 

in  terms  of  chosen  center  frequency  f  and  amplitude  and  phase 
modulations  A(t)  and  P(t),  respectively.  Although  f  and  P(t) 
are  not  unique,  the  argument  of  the  cosine  and  the  waveform  y(t) 
in  (A-7)  are  unique,  as  may  be  seen  by  the  first  equality  in 
(A-7).  All  of  these  relations  hold  for  time  t  e  T. 

If  the  fit  of  the  straight  line  a  +  pt  to  initial  unwrapped 
phase  P^(t)  over  interval  T  is  via  minimum  error  energy,  then  we 
find 


P  = 


U  V .  -  u 

Mo  1  o 

2  ' 

W2  -  "l 


n 


J  at  tn  , 

T 


vn  =  J  dt  tn  Pi(t)  .  (A-8) 

T 


There  is  no  need  to  explicitly  compute  a,  although  it  should  be 
included  in  the  error  energy  minimization  in  order  to  afford  a 
better  fit. 
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APPENDIX  B.  PROGRAM  FOR  FILTERED  COMPLEX  ENVELOPE  VIA  FFT 


The  program  listed  below  actually  computes  the  unfiltered 

complex  envelope  by  means  of  an  FFT.  In  order  to  convert  it  to 

one  which  will  compute  the  filtered  complex  envelope,  remove 

lines  220  -  320  and  replace  them  by  the  following  lines: 

220  Fl=40.  1  CUTOFF  FREQUENCY  <  Fo 

230  Hl  =  .  5*PI/F1 

240  Ml=M*Fl/Fo 

250  FOR  MS=-M1  TO  Ml  !  -FI  <  F  <  Fl 

260  J=Ms  MODULO  N  !  COLLAPSING 

270  F=Df *Ms  !  FREQUENCY  f 

280  CALL  Y ( F+Fo , A1 , Wo , Cp , Sp , Yr , Yi )  1SHIFTED  FREQUENCY  FUNCTION 

290  Cos=COS ( H1*F ) 

300  H=Cos*Cos  !  REAL  LOWPASS  HAWN  FILTER 

310  X( J)=X( J)+Yr*H 

320  Y ( J ) =Y ( J ) -Yi*H  1  CONJUGATE  INPUT  INTO  FFT 


10 

!  COMPLEX.  ENV 

ELOPE  VIA  SHIFTED  FREQUENCY  FUNCTION 

20 

fll  =1 

!  DAMPING  ALPHA 

30 

F  O  =  1  0  0 

i  CARRIER  FREQUENCY 

40 

Phi =-PI  -2 

!  PHASE 

50 

M  =  3  0  0  0 

!  NUMBER  OF  SAMPLES  FOR  F  <  0;  LINE  270 

60 

H  =  1  0  2  4 

!  SIZE  OF  FFT;  ZERO  SUBSCRIPT 

70 

REDIM  C  o i <  0 

:  N  •  ’  4 

> ,  X  <  0  :  N  -  1  >  ,  Y  <  0 :  N  - 1  ) 

30 

DIM  Cos'.  102 

4  >  ,  X 

<  4096) , Y (4096 ) 

90 

DOUBLE  M , N , 

N2,  J 

, Ms  !  INTEGERS,  NOT  DOUBLE  PRECISION 

100 

N2=N/2 

1  10 

fl  =  2 .  *  P  I  /  N 

120 

FOR  J=0  TO 

N/4 

130 

C  o  s  J  =  C  0  3 

A  *  J  > 

!  QUARTER-COSINE  TABLE  IN  Cos  CO 

140 

NEXT  J 

150 

C  p  =  C  0  3  <  P  h  i  > 

160 

3  p  =  3 1 N  ( P  h  i  > 

170 

Wo=2*PI*Fo 

130 

D  t"  =  F  o  /  M 

i  FREQUENCY  INCREMENT 

1  90 

Dr.  =  1  .  H * D  f 

>  TIME  INCREMENT  ON  COMPLEX  ENVELOPE 

200 

MAT  X='0.> 

2  1  0 

MAT  Y=  '  0  .  :■ 

220 

M  s  =  -  M 

230 

J  =  M s  MODULO 

N 

240 

CALL  Y  <  0 .  ,  A  1  ,  W o 

, Cp, Sp, Yr , Yi >  !  Y(0> 

250 

X  <  J  >  - .  5  *  Y  r 

260 

Y( J >=-. 5*Yi 

!  CONJUGATE  INPUT  TO  FFT 

270 

FOP  M s  =  - M  + 1 

TO 

M*10  i  NOTICE  UPPER  LIMIT  ON  FREQUENCY 

230 

J  =  M s  MODULO 

N 

!  COLLAPSING 

290 

F  =  D  f  *  M  s 

!  FREQUENCY  f 

300 

L  H  L  L  Y 'F  +  F o 

,  A  1  , 

Wo, Cp, Sp, Yr , Yi >  !  SHIFTED  FREQUENCY  FUNCTION 

3  1  0 

X  J  1  =  X  '  '  J  1  +  Y 

320 

■{<  J  >=Y  (  J  ■  -Y 

l 

!  CONJUGATE  INPUT  TO  FFT 

330 

NEXT  Ms 
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340 

MAT  X=X*<2. *Df > 

350 

MAT  Y=Y*<2.*Df> 

!  *  CONJUGATE 

360 

CALL  Fftl4<N,Cos< 

*>,X<*9, Y<*>>  !  y  <t)  OF  COMPLEX 

370 

G  I H  I  T 

!  -  ENVELOPE 

380 

PLOTTER  IS  "GRAPHICS" 

390 

GRAPHICS  ON 

400 

WINDOW  - N 2 , N 2 , - 8 , 

0  !  CENTER  PLOT  AT  TIME  t 

410 

LINE  TYPE  3 

420 

GRID  N/3, 1 

430 

LINE  TYPE  4 

440 

MOVE  0,0 

450 

Ts=Dt*N/4 

460 

DRAW  N/4, LGTCEXP< 

-  A 1  *  T  s  >  ) 

470 

PENUP 

430 

LINE  TYPE  1 

490 

FOR  M  s  =  -  N  2  TO  N2 

500 

J  =  M s  MODULO  N 

510 

X  =  X<.  J> 

520 

Y  =  Y  <  J  > 

530 

T  =  X  *  X  +  Y  *  Y 

'  MAGNITUDE  SQUARED  COMPLEX  ENVELOPE 

540 

IF  T : 0 .  THEN  570 

550 

PENUP 

560 

GOTO  580 

570 

PLOT  M=- ,  LGT '  T  5 

!  MAGNITUDE  OF  COMPLEX  ENVELOPE 

530 

NEXT  Ms 

590 

PENUP 

600 

PAUSE 

6  1  0 

GCLEAR 

620 

GRAPHICS  ON 

630 

WINDOW  - N 2 , N 2 , -PI 

,P  I 

640 

LINE  TYPE  3 

650 

GRID  NS, PI, 2 

660 

LINE  TYPE  1 

670 

FOR  M s  =  - N 2  TO  N2 

!  PLOT  COMPLEX  ENVELOPE  PHASE 

630 

J  =  M s  MODULO  N 

690 

PLOT  Ms, FNArg'XCJ 

>,-Y'J)>  i  CONJUGATE  THE  FFT  OUTPUT 

700 

NEXT  Ms 

7  1  0 

PENUP 

720 

PAUSE 

7  30 

GCLEAR 

"40 

GRAPHICS  ON 

750 

LINE  TYPE  3 

760 

GRID  NS, PI  2 

"70 

LINE  TYPE  1 

7  3  0 

FOP  Ms=-N2  TO  M2 

7  9  0 

J  =  M 1  MODULO  N 

3  0  0 

Ts=Ms-Dt 

'  TIME  r. 

3  1  0 

Cos =  COS  ■ .  Wo-Ts  ■ 

'  SHIFT  PHASE  THE 

3  2  0 

Si n= 3  IN' Wo-Ts  ■ 

'  COMPLEX  ENVELOFE  DY  Wo  Ts 

3  3  0 

X  =  X  (  J  1 

3  4  0 

Y  =  -  i'v  J 

'  CONJUGATE  THE  FFT  OUTPUT 

3  50 

PLOT  M s , F N A r  g 1  X - C 

o  s  -  Y  *  S  i  ri ,  X  *  S  i  n  +  Y  *  C  o  s  ) 

360 

next  Ms 

3  70 

PENUP 

3  3  0 

PAUSE 

390 

END 

3  0  0 

i 
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919  DEF  FHflrg < X ,  Y )  !  PRINCIPAL  ARG<Z!> 

920  IF  X=0 .  THEN  RETURN  ,5*PI*SGN<Y> 

930  ft=flTN('<VX> 

940  IF  X>0.  THEN  RETURN  FI 

950  IF  Y< 0 «  THEN  RETURN  fl-PI 

960  RETURN  A  +  PI 

970  FNEND 

980  ! 

990  SUB  Y<F,A1 ,Uo,Cp,Sp,Yr,Yi>  !  SPECTRAL  FUNCTION 
1000  W=2.*PI*F 

1010  T=W-Wo 

1020  D=AI*A1+T*T 

1030  R 1 = < Cp*A 1 +  Sp»T > ✓  D 

1O40  I  1  =  < Sp*A 1 -Cp*T ) /  D 

1050  T  =  W  +  W  o 

1060  D=A1*A1+T*T 

1070  R2= ( Cp*A 1 -Sp*T > /D 

1080  I  2  = ( - S  p  *  A  1  - C p  *  T  > / D 

1090  Yr  = . 5*  <  R 1 +R2  > 

1100  Yi=.5*<Il+I2> 

1110  SUBEND 

1120  i 

1130  SUB  F  f  t 1 4  <  DOUBLE  N , REAL  Co4<*>,X<*>, Y<*> >  !  N<  =2  1 4= 1 6334 ;  0  SUBS 
1140  DOUBLE  Log2r.,Nl,N2,N3,N4,  J,K  !  INTEGERS  <  2"31  =  2.147,483,648 

1150  DOUBLE  I  1  ,  12,  13, 14,  15,  16, 17, 13,  19, I  10, 1 1  1 ,  I  12.  I  13.  I  14, L' 0:  13. > 

1160  IF  N= 1  THEN  SUBEXIT 

1170  IF  N >2  THEN  1250 

1130  A  =  X  <  0 ) + X  <  1  > 

1190  X< 1 )=X(0>-X<  1  > 

1200  X  <  0  ‘  =  A 

1210  A  =  Y ( 0  >  +  V <  1  > 

1220  Y< 1 ) aY  <  0  > -Y  < 1 > 

1230  Y  <  0  >  =  A 

1240  SUBEXIT 

1250  A=L0G<N)/L0G<2. > 

1260  Log2n=A 

1270  IF  ABSC A-Log2n><  1 . E-8  THEN  1300 

1230  PRINT  "N  = " ; N ; "IS  NOT  A  POWER  OF  2;  DISALLOWED." 

1290  PAUSE 

1300  N 1  =  N  ■' 4 

1310  N  2  =  N 1  +  1 

1320  N  3  =  N  2  +  1 

1330  N  4  =  N  3  +  N 1 

1340  FOR  11  =  1  TO  L o g 2 n 

1 350  I  2  =  2  ■  < Log2n-  I  1  > 

1360  I  3  =  2+12 

1  370  I4  =  N,  13 

1380  FOR  15=1  TO  12 

1390  16=' 15-1 >* 14+1 

1400  IF  1 6 <  =  N 2  THEN  1440 

14  10  A 1  =  -  C  o  s  ( N  4  - 1 6  - 1  :> 

1420  A  2  =  - C  o  s ( 1 6 - N 1  -  1  > 

1430  GOTO  1460 

1440  Al=Cos< 16-1  ' 

1450  A 2  =  - C o s  <  N 3 -  I  6 -  1 : 

1460  FOR  17=0  TO  N-I3  STEP  13 

1470  13=17+15-1 

1430  19=13+12 
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1490 

ti=x< is> 

1500 

T  2  =  X  <  1 9  > 

1510 

T3=VU8> 

1520 

T4=Y< 19) 

1530 

fl3=T 1 -T2 

1540 

fl4=T3-T4 

1550 

X< I3>=T1+T2 

1560 

Y< I8)=T3+T4 

1570 

X< I9>=fll*A3-fl2*fl4 

1530 

Y  < I9>=fll*fl4+fl2*A3 

1590 

NEXT  17 

1600 

NEXT  15 

1610 

NEXT  11 

1620 

I l=Log2n+l 

1630 

FOR  12=1  TO  14 

1640 

l  <:  1 2  - 1  j  =  i 

1  650 

IF  I  2  L  o  q  2  n  THEN 

1670 

1660 

L(I2-1)®2A< 1 1-I2> 

1670 

NEXT  12 

1630 

K  =  0 

1  6  9  0 

FOR  11=1  TO  L(13) 

1  700 

FOR  12=11  TO  L'-.  12 

>  STEP  La  3) 

1710 

FOR  13=12  TO  Lai 

>  STEP  La  2) 

1  720 

FOR  14=13  TO  LU0 

>  STEP  L(ll) 

1  730 

FOR  15=14  TO  L<9> 

STEP  L (1 0 > 

1  740 

FOR  16=15  TO  L  3  > 

STEP  L < 9 > 

1  750 

FOR  17=16  TO  L  •:  7  > 

STEP  L < 8 > 

1  760 

FOR  13=17  TO  L 6 > 

STEP  L  <  7  > 

1  770 

FOR  19=13  TO  L ■ 5 > 

STEP  L  <.  6  > 

1730 

FOR  110=19  TO  L<4)  STEP  L<5) 

1  790 

FOR  111  =  110  TO  L  < 

3  ">  STEP  L  '  4  > 

1  3  0  0 

FOR  112=111  TO  L  < 

2  >  STEP  L  <  3  > 

1310 

FOR  113=112  TO  L ( 

1  >  STEP  L'12> 

1  320 

FOR  114=113  TO  L < 

0 )  STEP  L a  > 

1330 

J  = I 14-1 

1340 

IF  K >  J  THEN  1910 

1350 

fl  =  X  ■:  K  > 

1  Si  6  0 

X  <,  k  >  —  X  J 

13  70 

X  <  J  >  =  fl 

1330 

A  =  Y  ■'  K  > 

1  3  9  0 

Y <  K  1  =  Y J  > 

1  3  0  0 

Y  (  J  >  =  fl 

1  9  1  0 

K  =  K  +  1 

1  9  2  0 

NEXT  114 

1  9  3  0 

NEXT  113 

1  >4fl 

NEXT  112 

:-5  0 

NEXT  Ill 

1  >60 

NEXT  110 

1  970 

NEXT  19 

1  >30 

NEXT  13 

1  990 

NEXT  17 

2  0  0  0 

NEXT  16 

2010 

NEXT  15 

2O20 

NEXT  14 

20  3  O 

NEXT  13 

2  040 

NEXT  12 

'  0  5  0 

NEXT  11 

2  0  6  0 

SUBEND 
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APPENDIX  C.  CONVOLUTION  OF  TWO  WAVEFORMS 

Suppose  real  waveform  x(t)  excites  passband  filter  H(f)  with 
real  impulse  response  h(t).  Then,  the  output  is 

Y(f)  =  H(f)  X ( f )  ,  y ( t )  =  h(t)  ©  x(t)  .  (C-l) 

The  single-sided  output  spectrum  is 

Y+(f)  =  2  U ( f )  Y ( f )  =  2  U(f)  H(f)  X ( f )  =  \  H+(f)  X+(f)  .  (C-2) 

The  corresponding  output  analytic  waveform  is  exactly 

Y+(t)  =  \  h+(t)  ©  x+(t)  ,  ( C  —  3 ) 

which  is  just  (one-half  of)  the  convolution  of  the  individual 
analytic  waveforms. 

If  the  center  frequency  of  Y+(f)  is  fc  (see  appendix  A),  then 
the  spectrum  of  the  output  complex  envelope  is,  using  (C-2), 

Y(f)  =  Y+(f+fc)  =  \  »+(f+fc)  X+(f+fc)  =  \  H(f)  X(f)  ,  ( C-4 ) 

where  we  have  taken  the  same  center  frequency,  fc,  for  H+(f)  as 
well  as  X+(f).  This  relation  in  (C-4)  is  exact;  it  involves  no 
narrowband  approximations .  The  output  complex  envelope 
corresponding  to  (C-4)  is  then  exactly 

y(t)  =  \  h(t)  ©  x(t)  .  (C-5) 

That  is,  the  complex  envelope  of  the  convolution  of  any  two 
waveforms  is  equal  to  (one-half  of)  the  convolution  of  the  two 
individual  complex  envelopes,  irrespective  of  their  frequency 
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contents . 

Now  suppose  that  x(t)  is  given  in  terms  of  some  complex 
imposed  modulation  x^(t)  according  to 

x(t)  =  Re{x^(t)  exp( i2nf ct ) }  ,  ( C—  6 ) 

which  allows  for  amplitude-modulation  as  well  as  phase- 
modulation.  The  spectrum  of  x{t)  can  then  be  expressed  as 

X ( f )  =  \  [Xi(f-fc)  +  x*(-f-fc)]  .  ( C  —  7 ) 

Also,  suppose  that  filter  impulse  response  h(t)  is  expressible  in 
a  similar  form  according  to 

h(t)  =  Re{hi(t)  exp( i2nfct) }  ,  (C-8) 

with  corresponding  transfer  function 

H ( f )  =  \  [Hi(f-fc)  +  H*(-f-fc)]  .  (C-9) 

The  filter  output  spectrum  then  follows  from  (C-l),  (C-7), 
and  (C-9)  as 

Y<£>  -  j  tHi( f-fc)  Xi(£-fc)  +  H*(-£-£c)  X*(-£-£c)  + 

+  H*(-f-fc)  X.(f-fc)  +  Hi(f-fc)  X*(-f-fc)  ]  .  (C-10) 

By  inverse  Fourier  transforming  the  individual  terms,  the 
corresponding  waveform  to  (C-10)  is  found  to  be  exactly 

y ( t )  =  Re(exp( i2nfct)  [yg(t)  +  yb(t)])  ,  (C-ll) 

where 
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ya(t)  =  \  hi(t)  ©  xi(t)  ,  (C-12  ; 

and 

yb(t)  =  \  exP(~i4nfct)  3  ®  xi(t)  .  (C-13) 

Relation  (C-12)  states  that  component  y  ft)  of  output  y(t) 

a 

in  ( C— 11)  is  just  the  convolution  of  the  two  complex  imposed 
modulations  h^(t)  and  x^(t).  However,  (C-ll)  and  (C-13)  reveal 
that  there  is  an  additional  term  in  y(t),  which  requires  che 
convolution  of  a  relatively  high-frequency  component,  namely 
exp( -i4nf ct ) .  Since  this  latter  term,  y^t),  will  often  be 
Siaall  due  to  this  oscillatory  integrand,  we  may  neglect  it  in 
many  circumstances. 

A  good  way  of  assessing  the  importance  of  the  yb(t)  term  in 
(C-ll)  is  to  observe  that  it  is  due  to  the  second  line  of  the 
spectrum  in  (C-10);  the  first  line  in  (C-10)  corresponds  to 
y_(t).  Since  H-(f)  and  X-(f)  are  generally  lowpass  functions  of 
frequency,  the  function  H.(-f-f  )  in  (C-10)  is  centered  around 
f  =  -f  ,  while  the  X-(f-f  )  term  peaks  near  f  =  f  .  The 
separation  of  these  two  functions  is  approximately  2fc  on  the  f 
axis;  if  this  separation  is  somewhat  greater  than  the  bandwidths 
of  and  X^,  then  there  is  inconsequential  overlap  of  any  of  the 
frequency  components  in  the  second  line  of  (C-10).  This  leads  to 
a  small  value  for  y^(t)  for  all  t  and  we  can  neglect  its  effect 
relative  to  ya(t)  in  (C-ll). 

ct 
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